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1. SCOPE 


This document contains the "As-Built" Design Specifications for 
implementing Procedure 1 in EOD-LARSYS. Software for two new 
processors were written and routines for four existing processors 
were modified. 

EOD-LARSYS is operational on the Univac 1108 EXEC II version 
in Building 12. The batch oriented system is operated and 
maintained according to IDSD procedures. 

This document assumes the reader is familiar with both Procedure 1 
and the LARSYS system. 


2 . APPLICABLE DOCUMENTS 


• User Documentation EOD-LARSYS, Lockheed Electronics Company, 
Inc., HASD, Houston, Texas, November 1975 

• Job Order: 63-1557-1695 

• Requirements Document, Ref: 642-2240, dated February 7, 1977 

• TIRFi 77-0008 

• Design Specification for LARSYS Procedure 1, Lockheed 
Electronics Company, Inc., HASD, Houston, Texas, April 1977 



3. SYSTEM DESCRIPTION 


This document describes the EOD-LARSYS Procedure 1. In particular, 
this section is intended to outline the sequence of executing the 
various processors to accomplish the task of classifying a LACIE 
segment. These capabilities described are those of first priori- 
ty, and it is expected that other features will be added in the 
near future. The system is designed to provide RT&E with a means 
for experimenting with a new technique for classifying Landsat 
data for LACIE. 

The data is preprocessed by ERIPS. (ERIPS merges tapes received 
from Goddard Space Flight Center and creates a multitemporal/ 
multi-pass tape.) Using an ERIPS tape as input to the DOTDATA 
processor, a dot data file is output. The files contain both 
type 1 and type 2 dots. Type 1 dots are used both as starting 
vectors for the clustering processor (ISOCLS) and as labeling 
vectors for the labeling processor (LABEL) . Type 2 dots are used 
as a bias correction factor in computing the classification 
results output by the DISPLAY processor. 

After delineating DO/DU fields by card input, ISOCLS clusters the 
segment using the starting vectors from the dot data file to 
initialize the clustering process. Sun angle correction is pro- 
vided. An unconditional cluster map and a set of 'unlabeled' 
stats is output. 

The 'unlabeled' stats, cluster map, and dot data file are input 
to the LABEL processor. Using one of two procedures, k-nearest 
neighbor or all-of-a-kind, the stats are 'labeled'. A conditional 
or mixed cluster map may be output to be later displayed on the 
PMIS DAS or Image-100. 
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Using the 'labeled' stats, and allowing the user the capability 
of setting the inter-subclass weights by categories, the beft k 
of n channels may be selected by the feature selection processor 
SELECT. 

In CLASSIFY, each class in the 'labeled' stats may be assigned 
to a category by the analyst or by the system. The a-priori 
probability for each category may be computed using the cluster 
population from the stats. Using the k best channels, the sum-' 
of-densiti classifier assigns each pixel in the LACIE segment 
to a given subclass. 

In DISPLAY, the bias correction dots are used in computing the 
bias correction for. the classified area. Two performance tables 
relating to the bias correction dots are output. One table 
contains a category dot summary performance; the other, an indi- 
vidual dot summary performance. 

Given that the analyst is not satisfied with the classification 
results, he may relabel the dots and/or relabel the stats (by a 
labeling procedure or by card input) and restart the process at 
any given point. 

3 . 1 HARDWARE DESCRIPTION 


N/A 


3.2 SOFTWARE DESCRIPTION 


3.2.1 SOFTWARE COMPONENT NO. 1 (DOTDATA) 

In implementing Procedure 1 , the means for allowing the user the 
caped^ility to l 2 d)el certain MSS data points (known as dots) was 
added to the system by the processor DOTDATA. The main function 
of this processor is to output a file containing the dots of 
interest. This file is an interface for three processors, ISOCLS, 
LABEL, and DISPLAY. 

The dots are defined by field cards. By user control, any subset 
of the 209 possible grid points may be selected. Optionally, 
the dots may be labeled at the initial creation of the file or 
the dots may be labeled in the LABEL processor. All the cate- 
gories of interest do nol. have to be defined in this processor. 

If a dot is to be retyped, tihe file must be recreated by executing 
this processor again. 

The files contain both type 1 and type 2 dots. Type 1 dots 
(starting and labeling dots) are written on one file; type 2 
dots (bias correction dots), on a second file. 

By an OPTION control card, the user may request that the spatial 
and spectral information relating to each dot on the file be 
printed on the line printer. 

3. 2. 1.1 Linkages 

The processor DOTDATA uses the FORTRAN-V compiler, Uni vac soft- 
ware system routines, EOD-LARSYS utility routines, and common 
blocks INFORM, GLOBAL and DOTVEC, and the following processor 
subprograms: DOTDAT, DOTS, FLDTYP, and SET13- 


See Appendix A for description of subprograms . 


3. 2. 1.2 Interfaces 


N/A 


3. 2. 1.3 Inputs 

Requires an MSS data tape. See section 3.2 of User Documentation 
EOD-LARSYS. Requires three types of card input. 

• Processor Card 


Keyword Function 

{Col. 1) 

$DOTDATA Loads into the system 

all the routines needed 
for executing this pro- 
cessor. 


• New Control Cards 


Keyword 
(Col. 1) 
CHANNEL 


DATAFILE 


DOTFIL 


Parameter 

(Col. 11-72) 

DATA=n^,n2, • • • ,n^Q 
(Default: None) 

UNIT=n,FILE=m 
(Default: n=3 
m=l) 


OUTPUT/UN I T=n , FILE=m 
(Default: n=8 
m=l) 


Function 

Integer numbers separated 
by commas referring to 
the channels on the MSS 
data tape. 

n is the FORTRAN unit 
number assigned to the 
MSS data tape, m is 
the file number of the 
data to process. 

n is the FORTRAN unit 
number assigned to the 
dot data file output by 
this processor, m is the 
number of the file to 
output. 




Keyword 

Pareuneter 

Function 

OPTION 

PRINT 

(Default: no line 
printer output) 

Initiates the printing 
of the dot data file 

information. 

HEDl 

Any 60 characters 
(Default: Lyndon 
B. Johnson Space 
Center) 

First line of the 
heading on line 
printer output. 

HED2 

Any 60 characters 
(Default: Houston, 
Texas) 

Second line of the 
heading on line 
printer output. 

DATE 

Any 12 characters 
(Default: Present 
date) 

Date in the heading 
on line printer out- 
put. 

COMMENT 

Any 60 characters 
(Default: None) 

Comment printed with 
heading on line printer 
output. 

★END* 


Indicates the end of 
the control card inputs. 

$END* 


Indicates the end of 
all the card inputs. 


Field Cards 

By the field card, the user defines the grid ^joints to 
extract from the MSS data tape. The definition and order of 
the field card(s) determines the position of the dots on the 
dot data file, DOTFIL. The analyst will need to know the 
position of the dots for defining starting vectors in ISOCLS 
and for labeling/relabeling the dots in LABEL. 


At the time of defining the fields, the type for each dot is 
defined by the TYPE card. By option, the analyst may label 
each dot by a CLASSNAME card. If this card is omitted, the 
unlabeled dots must be labeled by the control card DOTLABEL 
or excluded from the set by the control card EXCLUDE in the 
labelling processor, LABEL. 

An example is given to illustrate a field data set expected by 
this processor: 

*END* 

TYPE 1 

CLASSNAME WHEAT (optional) 

LABELl (10,10) , (10,10) , (196,10) 

LABEL2 (10,10) , (10,20) , (196,20) 

CLASSNAME NONWHT (optional) 

LABELS (10,10) , (1C, 50) , (100,50) 

TYPE 2 

CLASSNAME WHEAT (optional) 

BIASl (10,10) , (10,40) , (196,40) 

CLASSNAME NONWHT (optional) 

BIAS2 C ^,10) , (10,70) , (196,70) 

$END* 

Two files are written. File 1 contains 38 WHEAT dots, followed 
by 10 NONWHT dots. All of which are type 1 dots. File 2 contains 
19 WHEAT dots followed by 19 NONWHT dots . All of which are 
type 2 dots . 

If the CLASSNAI'lE cards were omitted, file 1 would contain 48 
unlabeled, type 1 dots. File 2 would contain 19 unlabeled 
type 2 dots . 

In both cases, the reference number for the dots in file 1 
defined by LABELl field card is 1 thru 19; LABEL2 field card 
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is 20 thru 38; and LABELS field card is 39 thru 48. The refer- 
ence number for the dots in file 2 defined by BIASl field card 
is 1 thru 19 and by BIAS2 field card is 20 thru 38. 

3. 2. 1.4 Outputs 

DOTFIL, a multi-file unformatted FORTRAN written tape, is output 
as an interface to the processors ISOCLS, LABEL, and DISPLAY. 

The default unit for DOTFIL is logical unit F (FORTRAN unit 8) . 
See Appendix H for format of the tape. 

The line printer output is: 

• Summary of field information 

• Spatial and spectral information for each dot (optional) 

• Summary of user options 

3. 2. 1.5 Storage Requirements 
Requires less than 53K of core. 

3 . 2 . 1 . 6 Description 

The main function of DOTDATA is to create the dot data file, 
DOTFIL. Optionally, the spatial and spectral information of 
each dot in the DOTFIL may bo printed on the line printer. 

DOTDAT is the driver routine. From this routine, the setup 
routine SET13 is called. SET13 reads and analyzes the control 
cards and prints a summary of the user selected options. 

DOTS is the coordinator and controller of the routines needed 
to create and output the file, DOTFIL. 

From DOTS, FLDTYP is called to initiate the reading of the field 
cards. In returning from FLDTYP, the sample number, lino number, 


type, category (optional) , and data vector are stored in the 
array DATA. DATA is passed to the routine DOTFIL; this routine 
outputs the dot data file. 



3. 2. 1.7 Flowchart 






3.2.2 SOFTWARE COMPONENT NO. 2 (LABEL) 

To aid the analyst in supervising the labeling of the statistics 
obtained from the clustering processor ISOCLS, a now technique 
for labeling the statistics was implemontv .. 

Two procedures for labeling the statistics is provided. 

Optionally, the analyst may select either the k- nearest neighbor 
procedure or the all-of-a-kind procedure. 

The labels in the dot data file, DOTFIL, and/or previously 
labeled statistics file, SAVTAP, may be relabeled by control 
card input and the updated file(s) output to tape. Optionally, 

1) a conditional or mixed cluster map may be output to taoe, 

2) an unconditional cluster map may be output in the format 
acceptable by the DISPLAY processor, 3) the spatial and spectral 
information of the relabeled DOTFIL may be output to the line 
printer, and 4) the statistics of the SAVTAP file may be output 
to the line printer. 

3. 2. 2.1 Linkages 

The processor LABEL uses the Fortran-V compiler, Univac software 
system routines, EOD-LARSYS utility routines, common blocks INFORM, 
GLOBAL, and LABS, and the following processor subprograms: LABEL, 

KNEAR, ASCEND, SET14 , CRDSCN, FILERD, STOMAP, LABLR, DSPTAP, 

MIXMAP, CNDMAP, NAMSTA, LABDOT, LABMAN, REDDER, MANORD, DOTDST, 
SUNFCS, ALLKIN, CLSMAP, CLRKEY, and MAPHND. 

See appendix B for description of subprograms. 

3. 2. 2. 2 Interfaces 

Optionally, the three files listed below may be input: 

• Statistics file (SAVTAP) from ISOCLS or STAT processor 

• Cluster map file (MAPFIL) from ISOCLS processor 

• oot data file (DOTFIL) from DOTDATA processor 



For format of the files see section 4.1 and 5.1 of the User 
Documentation EOD-LARSYS and appendix H of this document. 


3. 2 . 2. 3 Inputs 

Requires 3 types of card input. 

• Processor card 

Keyword Parameter Function 

(coi. 1) 

$LABEL Loads into the system all the 

routines needed for executing 
this processor 

• New control cards 


Keyword 
(col. 1) 
CHANNEL 


Parameter Function 

(col. 11-72) 


STAT=nj^,n2,-“,n3p. 

DATA^rn^^ , m2 » * • * » m^Q 

(Default; 

n^ = {all channels 
on SAVTAP 
file} 

m^ = (all channels 
on DOTFIL 
file} 


n^ and m. are integer numbers 
separated by commas referring 
to the channels on the SAVTAP 
file and the DOTFIL file, 
respectively. 


DOTFILE INPUT/FILE=m,UNlT=n Defines the unit and file for 

the input dot data file, DOTFIL. 
n is the Fortran unit number 
assigned to the input DOTFIL 
file, m is the number of the 
input file to process. 



Function 


Keyword 

DOTFILE 


STATFILE 


STATFILE 


MAPTAP 


Parameter 

OUTPUT/FILE=m, UNIT=n 
(Default: NONE) 


INPUT/FILE=m, UNIT=n 
(Default: NONE) 


OUTPUT/FILE=m, UNIT=n 
(Default: NONE) 


OUTPUT/FILE=m , UNIT=n 
(Default: no DISPLAY 

interface tape will be 
output. ) 


Defines the unit and 
file to which the re- 
labeled DOTFIL file 
is output. 

n is the Fortran unit 
number assigned to the 
output DOTFIL file, 
m is the nvimber of the 
file to output. 

Defines the unit and 
file for the input 
SAVTAP file, 
n is the Fortran unit 
number assigned to 
the input SAVTAP file, 
m is the number of the 
file to process. 

Defines the unit and 
file to wnich the labeled/ 
relabeled SAVTAP file 
is output. 

n is the Fortran unit 
number assigned to 
the output SAVTAP file, 
m is the number cf 
the file to output. 

Defines the unit and 
file to which the 
unconditional cluster 
map, MAPTAP, is output, 
n the Fortran unit 
number assigned to the 
output MAPTAP file. 




Keyword 


Pariuneter 


Function 


MAPFILE 


MAPFILE 


DOTLABEL 


INPUT/FILE=m,UNIT=16 
(Default; m*l) 


OUTPUT/FILE=m, UNIT=n 
(Default: n=16 

m=l 


category name, 

”l'”2' * ' * '^250 
(Default; NONE) 


If executing back to back 
with DISPLAY, n must be 
assigned to unit 2. m is 
the number of the file to 
output. 

Defines the unit and file 
for the input unconditional 
cluster map, MAPFIL, that 
is output by ISOCLS during 
the clustering process. 
n=16 is the Fortran unit 
number assigned to the 
input MAPFIL file. (n 
must be assigned to unit 16 
if executing back to back 
with ISOCLS) . m is the 
number of the file to 
process. 

Defines the unit and file 
to which the conditional 
and/or mixed cluster map is 
output, n is the Fortran 
unit number assigned to 
the output MAPFIL file, 
m is the number of the 
file to output.^ 

The DOTFIL is labeled or 
relabel by this card. Cate- 
gory name is the label the 
analyst is assigning to the 


If both type of maps are output, the conditional map is out- 
put on file m; the mixed map, on file m+1. 



Keyword 


Function 


STATLABEL 


DISTANCE 


Parameter 


class name, n. , 

1 

n2,“-,n25o 

(Default: NONE) 


LI 

(Default: Li 

distance) 


dots n j , j“l,**',i. The 
category name may be 
composed of a maxi- 
mum of 6 characters, 
n^ are integer numbers 
separated by commas 
referring to the posi- 
tion of the dot on 
the DOTFIL file. 

The SAVTAP file may 
be manually relabeled 
by this card, 
n j (3=1 , 2 ... X ) are the 
number of the sub- 
classes on the SAVTAP 
that are to be re- 
grouped into another 
class. Class name is 
the name of the class 
to which the subclasses 
nj are to be reassigned. 
Class name must match 
a name on the SAVTAP 
file. 

The L]_ distance is 
used in computing the 
distance between the 
means of the cluster 
and the labeling 
dots. 


/y 


Keyword 

Parameter 


Function 

DISTANCE 

L2 


The distance is 


(Default; 

distance) 

^1 

used in computing 
the distance between 
the means of the 
cluster and the 
labeling dots. 

OPTION 

COND 


A conditional cluster 


(Default; 

NONE) 

map i s output . 

OPTION 

MIXED 


A mixed cluster 


(Default; 

NONE) 

map is output. 

THRESHOLD 

T 


T is the threshold 


(Default: 

T=25.0) 

parameter used in 
creating the condi- 
tional cluster map. 

T is a floating point 
number. 

NEAREST 

K 


K is the number of 


(Default: 

K=l) 

dots to be used in 
the k-nearest neights 
bor procedure. K is 
an integer number. 
K<250 . 

PROCED 

NAME 

(Default: 


NAME is an alpha word. 
NAME=K-NEAREST (Use the 


N=K-NEAREST 

k-nearest neighbor 


procedure. ) 

NAME=ALL (Use the all- 
of-a-kind procedure.) 


Keyword 


Function 


MODULE 


EXCLUDE 


SUNANG 


SUNANG 


Parameter 


(blank) 


* * ’ '^250 

(Default: all dots 

on the DOTFILare 
used. ) 


(Default: no sun 

angle correction is 
be applied.) 


PILE 

(Default: No sun 

angle correction 
is applied.) 


NAME-MANUAL (Use the manual 
procedure of relabeling 
the DOTFIL or SAVTAP file. 

Initiates the input of 
the module STAT card 
deck. The deck must 
immediately follow this 
card. 

n^^ are integer numbers 
referring to the dots on 
the DOTFIL that are to 
be excluded in all calcu- 
lations (i.e., dots 
within a DO/DU area) . 

m^ are integer sun angle 
numbers used in computing 
the Lj^ or L 2 distances. 

A sun angle must be input 
for each acquisition of 
interest. An acquisition 
is assumed to be a 4 chan- 
nel pass. 

Example: If the distance 

is computed using 16 chan- 
nels, 4 sun angles (m^, 
m 2 , m^f m^) must be input. 

Sun angles will be 
extracted from the 
DOTFIL. 
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Keyword 

Pareuneter 

Function 

OPTION 

STATS 

(Default: Stats not 

printed) 

Means and covariances for 
labeled or relabeled 
statistics on SAVTAP file 
are printed. 

OPTION 

DOTS 

(Default: Relabeled 

DOTFIL is not printed) 

Spatial and spectral 
information of relabeled 
DOTFIL is printed. 

HEDl 

Any 60 characters 
(Default: Lyndon 

B. Johnson Space 
Center) 

First line of the heading 
on line printer output. 

HED2 

Any 60 characters 
(Default: Houston 

(Texas) 

Second line of the heading 
on line printer output. 

DATE 

Any 12 characters 
(Default: Present 

date) 

Date in the heading on 
line printer output. 

COMMENT 

♦END* 

Any 60 characters 
(Default: NONE) 

Comment printed with 
heading on line printer 
output . 

Indicates the end of the 
control card inputs. 




Keyword 

$END* 


Pariuneter 


Function 


Indicates the end of 
all the card inputs. 


• FIELD CARDS 

If 1) the all-of-a-kind procedure is selected, 2) a conditional 
or mixed cluster map is output, or 3) a DISPLAY interface tape 
is output, a MAPPIL tape must be input, and a field card 
defining the area of the unconditional cluster map (input 
MAPFIL) must also be input. The vertices must reflect the 
sample number and line number of the MSS data tape used to 
create the unconditional cluster map. That is, the field card 
must be the identical to the field card input to ISOCLS. 

If 1) a procedure is not to be executed, 2) a MAPFIL file not 
output, or 3) a MAPTAP file is not output, a field card is not 
input. 

Example of the data set; 

• MAPFIL tape is being input 
*END* 

FIELD (1,1), (1,1), (196,1), (196,117), (1,117) 

$END* 

• MAPFIL tape is not being input 
*END* 

$END* 

3. 2. 2. 4 Outputs 

Optionally the following files are output: 

• Labeled statistic file (SAVTAP) by using a procedure (See 
section 4.1 of User Documentation EOD-LARSYS for format of 
file) 


• Labeled statistic file (SAVTAP) by control card input 

• Relabeled dot data file (DOTTIL) (See Appendix G for for- 
mat of file) 

• A conditional cluster map (See section 5.1 of User Docu- 
mentation EOD-IiARSYS for format of file) 

• A mixed cluster map. 

• An unconditional cluster map in the format acceptable by the 
DISPLAY processor (See Appendix C of User Documentation 
EOD-LARSYS for format of file) 

i. xine printer svimmary of the following is output: 

• Summary of selected options 

• Table of L^ or L 2 distances 

• Summary of the labeling dots within a cluster for the 
all-of-a-kind procedure 

• Summary of the labeling dots for the k nearest dots to a 
cluster for the k-nearest neighbor procedure 

• Spatial and spectral information of relabeled DOTFIL 

• Means and covariances of labeled or relabeled statistics. 

3. 2. 2. 5 Storage Requirements 

Requires about 53k of core. 

3. 2. 2. 6 Description 

A distance table containing the or L 2 distance of all the 

type 1 dots and the mean of the clusters is computed. 



^'0 

'' 


where ; 


X. = i element of the dot vector 

X 

element of the mean vector 
n = number of channels 

Using the table as input to the k-nearest neighbor procedure or 
all-of-a-kind procedure, the statistics generated during cluster- 
ing is labeled. 

For the k-nearest neighbor procedure, all the labels of the 
k-nearest labeling dots to a given cluster are polled. The 
label with the majority of the dots labels the cluster. If a 
tie occurs, then k-1 dots are considered. 

For the all-of-a-kind procedure, all of the labeling dots within 
a cluster are polled. If all the dots are of one category, the 
cluster labels that category. If the cluster contains dots for 
more than one category, the label with the majority of the dots 
labels the cluster. If there are no labeling dots within a 
cluster, the labeling defaults to the k-nearest neighbor proce- 
dure. 

Optionally, a conditional olustez map may be output. A cluster 
is tagged as conditional if the distance between the nearest 
identically labeled labeling dot and mean of the cluster is 
greater than the analyst input threshold value t. 

Optionally, a mixed cluster map may be output. A cluster is 
tagged as mixed if the labeling dots within a cluster are of more 
than one category . 

Optionally, an unconditional cluster map may be output in the 
format acceptable to the DISPLAY processor. Information used in 




the thresholding procedure in DISPLAY is dummied. If thi ^sholding 
of the clustered data is desired, it can be performed by exercising 
the conditional map option in this processor. 
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3.2.3 SOFTWARE COMPONENT NO. 3 (ISOCLS) 

Several additions were made to the ISOCLS (clustering) processor 
of the LARSYS system in support of Procedure 1 requirements. 

These additions are described in a general sense as follows; 

a. On option starting dots (pixels) from the dot data file 
DOTPIL may be used to begin clustering. 

b. The analyst may identify "designated other" and "designated 
unidentifiable" pixel sets (fields) by field card input. 

The pixels in these fields are not included as inputs to the 
clustering algorithm. They are assigned special cluster 
numbers and mean vectors for display purposes . 

c. On option, using a sun angle correction table, the pixel 
radiance values may be modified. [The correction table is 
built in.] The radiance value correction applies only for 
clustering purposes. 

The user may input the sun angles by cards or request that 
these angles be extracted from the header record of a 
universal formatted MSS data tape (ERIPS unload tape) . 

3. 2. 3.1 Linkages 

The processor ISOCLS uses the FORTRAN-V compiler, Univac software 
system routines, EOD- LARSYS utility routines and common blocks 
PASS, GLOBAL and ISOLNK. 

See Appendix C for description of the subprogram SUNFAC and modi- 
fications to the following subprograms: PSPLIT, SUTUP7, RDDATA, 

and ISODAT. 

3 . 2 . 3 . 2 Interfaces 

Files generated by other processors and used as input to ISOCLS 
are: 






a. On option, a statistics file (SAVTAP) to provide starting 
cluster mean vectors. 

b. On option, input initial cluster "centers" by cards. 

c. On option, startiiig dots from DOTFIL to initialize the 
cluster processing. 

Format descriptions of those files are included in the User Docu- 
mentation EOD-LARSYS in sections 4.1 and 3. 1.4. 3, and Appendix H 
in this document, respectively. 

3. 2. 3. 3 Inputs 

Require an MSS data tape and card input. See section 3.2 
and 9.5, respectively, of User Documentation EOD-LARSYS. 

Added control cards and field cards are note below: 


• New Control Cards 


Keyword 


Parameters 


Function 


DOTFILE 


INPUT/UNIT=n,FILE=m Defines the FORTRAN unit 

(Default: Self- number n and file number m 

initializing starting) of the dot data file DOTFIL 

containing the starting 
vectors . 


DOTS 


SUNANG 


Hi '^ 2' ■ • ' '^60 

(Default: Dots will 

not be used for start- 
ing vectors.) 


n^ are integer numbers 
separated by a comma speci- 
fying the dots to be used 
as starting vectors. 


TAPE 

(Default: No sun 

angle correction 
applied . ) 


Sun angles a]:e extracted 
from the ERIPS unload MSS 
data tape. 


siy 


Keyword 


Function 


Parameters 

SUNANG n^^ , n2 » . . . n j 

n^ are integer 
nvimbers , j 8 . 
(Default: No 

sun angle correc 
tion applied. ) 


n^ are the sun angles to 
be used in computing the 
sun angle corrections for 
use in the clustering 
algorithm. A svin angle 
must be input for each set 
of 4 channels input on the 
CHANNEL control card. 


• Field Cards 


ISOCLS recognizes DO/DU fields. All the DO/DU field cards 
(for all classes) must be input before the fields to be 
clustered. These fields must immediately follow the *END* 
card. The CLASSNAME card follows the last DO/DU field card. 


Exeunple: 


• If DO/DU fields are being defined 
*END* 


DESIGNATED 

OTHER 

DESIGNATED 

UNIDEN 

CLASSNAME 

WHTl 

$END* 


OTHER 

(1.1) , (1,1) , (40,1) , (40,20), (1,20) 
UNIDENTIFIABLE 

(1.1) , (5,7) , (8,7) , (8,10) , (5,10) 

WHEAT 

(1.1) , (1,1) , (196,1) , (196,117) , (1,117) 


• If no DO/DU fields are being defined 
*END* 

CLASSNAME WHEAT 

WHl (1,1) , (1,1) , (196,1) , (196,117) , (1,117) 

$END* 



punched card deck 


3.2. 3. 4 Outputs 

a. Statistics file (SAVTAP) and on option, 

b. Cluster map tape (MAPFIL) , on option 

c. Printout of cluster results, consisting of the following 
data items by class: 

Cluster numbers and symbols 
Cluster mean vectors (by channel) 

Cluster standard deviations by channel 
Inter-cluster distances 

Number of pixels per cluster 
Nvunber of clusters 

Cluster map by field for each class 

3. 2. 3. 5 Storage Requirements 

Require 17449 words of core plus 19345 words of core for the 

resident programs. 

3. 2. 3. 6 Description 

The implementation of the additions to ISOCLS are generally 

described. For further details, please refer to the subprogram 

documentation which is included. 

a. By expanding the INVEC data vector of the control card is 
accepted. Specification of unit and file number on this 
card triggers a call to routine RDDOTS to load the 
starting vectors. These vectors are stored internally for 
use by ISODAT. The unit and file number is stored in an 
extension to labeled common block GLOBAL. 

b. Pixel fields to be described as "other" or "unidentifiable" 
are input by analyst-furnished classname and field definition 
cards to be read by RDDATA. The classname cards are punched 
as : 



DESIGNATE OTHER or 
DESIGNATE UNIDENTIFIABLE. 

The DO/DU pixels are flagged and assigned special cluster 
numbers. They are not included as input to the clustering 
algorithm. This requires a modification to RDDATA and ISODAT. 

Field definition information is stored internally for use by 
ISODAT. A trigger indicating that DO/DU fields have been 
used is set and stored in an extension of the labeled common 
block PASS. 

c. To handle sun angle input, SETUP? is modified (by adding 
•SUNANG' to the INVEC data vector) to read a sun angle 
control card. This card either contains in the parameter 
field of the card, the sun angles for each pass or the word 
TAPE. TAPE triggers the unpacking of sun angles from the 
MSS data tape. If the SUNANG control card is not present, a 
default to no sun angle correction occurs. The trigger 
indicating that the sun angles have been unpacked appears 

in labeled common ISOLNK. 

d. The sun angle correction table is stored as DATA in subpro- 
gram SUNFAC. The correction is computed by 1) looking up 
the gain versus sun angle for each channel, and 2) multi- 
plying the pixel radiance value for the channel by the gain. 
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3.2.4 SOFTWARE COMPONENT NO. 4 (SELECT) 

The SELECT processor is modified to optionally provide an auto- 
matic assignment of interclass weights. The weights for class^, 

class. = 1.0 for i 5 ^ j and the weights for class., class. = 0,0 
D ID 

for i ® j. The breakdown of clasr?^, clasSj pairs into the 

correct set of intersubclass pairs is provided by the processor. 

Information concerning the class-subclass association is extracted 

from the input statistics file, SAVTAP. 

The modifications to the SELECT processor involve two sxibprograms , 
SETUP4 and WGTCHK. These modifications are made to enable use of 
the processor in the evaluation of the LACIE Procedure 1 crop 
classification technique. For application of the SELECT processor 
to LACIE Procedure 1, the user must consider classes equivalent 
to categories. 

3. 2. 4.1 Linkages 

The SELECT processor uses the Univac 1108 (EXEC2) operating 
system and associated system routines, the FORTRAN V compiler, 
EOD-LARSYS utility routines, and EOD-LARSYS common blocks FSL, 
GLOBAL, and INFORM. 

See Appendix D for the detailed description of modifications to 
the subprograms SETUP 4 and WGTCHK. 

3. 2. 4. 2 Interfaces 

The SELECT processor interfaces with the LABEL processor via the 
labeled statistics file, SAVTAP, which is output by the LABEL 
processor. The statistics file is also the interface between 
SELECT and the STAT or ISOCLS processors. 

3 . 2. 4 . 3 Inputs 

The modifications to the SELECT processor include an additional 


input option on the OPTION control card. The new input option 
is required to initiate the automatic (processor provided) assign 


ment of interclass weights. See 
of the new option. 

• Revised Control Card 

Keyword Par^lmeter 

(col. 1) (col. 13-72) 

OPTION CLSWT 

(Default: The 
weights are 
assigned to 
intersubclass 
pairs. ) 


below for the format and function 


Function 

The processor determines the 
class-subclass correspondence 
(after any grouping of sub- 
classes to form new subclasses 
if the GROUP control card 
is used) and assigns a weights 
1.0 to the subclass pairs 
associated with all interclass 
pairs. Intraclass subclass 
pairs are given a weight=0.0 

NOTE(l): The WEIGHTS control 

card remains available to allow 
the user to set weights for 
specific subclass pairs. If 
used, the input subclass pair 
weights override the processor- 
set subclass pair weights. 

NOTE (2): The "WEIGHTS OTHERS" 

capability is not available 
when this option is exercised. 
If input, it is ignored by the 
processor. 


33 


3. 2. 4. 4 Outputs 

The message for automatic assignment of interclass weights is 
added to the line printer output of user requests. There is no 
change to file output by the SELECT processor. 

3. 2. 4. 5 Storage Requirements 

The current storage allocation (overlay) map for EOD-LARSYS 
results in a utilization of 7696 storage locations in the 
Instruction Bank, and 15309 locations in the Data Bank, plus 
19345 locations for the resident programs. 

3. 2. 4. 6 Description 

Initialized via the OPTION CLSWT control card, SETUP4 sets an 
internal flag (WTKEY=1) indicating that subclass pair weights 
are to be program-determined on an interclass basis. From the 
class (or category) and subclass information from the input 
statistics file, SAVTAP, the subclass-to-class association will 
13 determined. In determining the class-subclass association, 
any new subclasses created by the user via the GROUP control 
card is re*.cognized. 

Using the existing storage location, SUBRAY, SETUP4 builds a 
table of subclass pair weights. The weight table is initialized 
by a newly added internal subroutine, INTWGT. All subclass 
pairs which represent an interclass pair are given a weight=-1.0. 
All other subclass pairs (i.e., intraclass subclass apris) are 
given a weight=0.0. 

As usual, the user may input any subclass pair weights via the 
WEIGHTS control card. Any subclass pair weights input via the 
WEIGHTS card overrides the processor-set weights for those 
subclass pairs. This function is provided by the existing sub- 
program, WGTCHK, as usual. 
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Concluded - Select Processor Flowchart 



3.2.5 SOFTWARE COMPONENT NO. 5 (CLASSIFY) 

The classification processor, CLASSIFY, is modified as follows; 

1. To allow an option for obtaining subclass a priori values 
using subclass population data from the input file, SAVTAP. 

2. To allow the system tc assign the category names using the 
class names from the irput statistics file, SAVTAP, as the 
assigned category names. This modification is to adopt the 
processor to LACIE Procedure 1 techniques. 

Both options are an addition to the usual capability of analyst- 
input of a priori probability values at the subclass, class, or 
category level via the APRIORI control card and of category name 
input via the CATEGORY control card, 

3. 2. 5.1 Linkages 

The CLASSIFY processor uses the Univac 1108-EXEC 2 operating 
system and associated system routines, the FORTRAN V compiler, 
and EOD-LARSYS utility routines and common blocks. See Appen- 
dix E for the description of the modified subprograms. 

3. 2. 5. 2 Interfaces 

There are no modifications to CLASSIFY interfaces. As usual, the 
CLASSIFY processor obtains the class names, subclass statistics 
(mean vectors and covariance matrices) and related subclass data 
from either an input statistics file, SAVTAP, or the Module STAT 
card deck which are created by either the STAT processor or the 
ISOCLS processor. The CLASSIFY processor also interfaces with 
the SEnECT processor to obtain a transformation matrix (the B- 
matrix) via either a file, BMFILE, or a 3-MATRIX card deck which 
are both created by the SELECT processor. 
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The input statistics file, SAVTAP, is either assigned to UNIVAC 
unit A or to the users choice of available units. The input 
transformation matrix file, BMFILE, is always assigned to UNIVAC 
unit H. 


3 . 2 . 5 . 3 Inputs 

• MSS Data 

The MSS (multi-spectral scanner) data to be classified is 
input to the CLASSIFY processor via the input file, DATAPE, 
as usual. The file must be in either of two specific formats, 
LARSYS II or Universal, and is either assigned to Univac unit 
C, or to the user's choice of units as designated on the 
DATAFILE control card. 

• Revised Processor Control Card 

Changes to the control card input for the CLASSIFY processor 
to initiate the modifications described in Section 3.2.5 are 
on the APRIORI control card and CATEGORY control card. Either 
the modified or standard cards may be input, depending on the 


options 

desired. 


Keyword 

Parameter 

Function 

(col. 1) 

(col. 11-72) 


CATEGORY 

FILE 

(Default: No 

categories are 
defined and the 
standard classifier 
will be applied) 

Initiates the assigning 
of the category names 
using the class names 
from the input statistics 
file, SAVTAP, and invokes 
the category classifier. 


Function 


Keyword 

APRIORI 


Parameter 

FILE 

(Default: Subclass 

a priori is not 
be computed from 
the statistics file, 
SAVTAP. 


The subclass apriori 
probability values are 
computed using subclass 
or cluster point popula- 
tions from the statistics 
file, SAVTAP. 


• Field definition cards 

There are no modifications to the field definition input to 
CLASSIFY. As usual, the area(s) on the input data tape to be 
classified by the CLASSIFY processor are defined to the 
processor on input "field definition" card(s). The area(s) 
to be classified are defined in terms of sample, line coordi- 
nates of each vertex of the field, up to a maximum of ten (10) 
vertices for an irregularly shaped field. Also contained on 
the field definition card(s) are the incrementation of samples 
and/or lines to be performed in reading the data for the 
defined field. 


The format of the "field definition" card is provided in 
Section 3.1.3 of the user documentation EOD-LARSYS docu- 
ment no. LEC-3984. 


3 . 2 . 5 . 4 Output 

There are no modifications to the output provided by CLASSIFY. 

As usual, the classification processor, CLASSIFY, outputs the 
classification results on a file, MAPTAP, assigned to Univac 
Unit B. The primary results output on the MAPTAP file are the 
subclass identification of each pixel of the area classified, and 
the value of the probability density function for each pixel of 
the area classified. 




The classification processor outputs via the line printer a 
classification map, which illustrates the classified identity 
of each pixel by symbols, in correct spatial relationship according 
to seunple, line coordinates of each pixel. Also output to the 
line printer are the class/subclass training field coordinates, 
the class/subclass statistics (mean vectors and covariance 
matrices) and a listing of the subclasses considered for classi- 
fication along with apriori probability values for each subclass. 

3. 2. 5. 5 Storage Used 

With the present EOD-LARSYS overlay (MAP) structure, the CLASSIFY 
processor requires 4058 locations in the Instruction Bank, 24 55*^ 
locations in the Data Bank, for a total of 28617 locations, plus 
19345 locations for resident programs. 

3. 2. 5. 6 Description 

To accomplish the additional option described in Section 3.2.5, 
the modifications to the CLASSIFY processor involve two sub- 
programs, REDIF2 and SETUP2. 

Subprogram REDIF2 is modified to read and decode the revised 
APRIORI control card and CATEGORY control card described in 
Section 3.2.5. 3. If the word FILE is detected in the parameter 
field (card columns 11 through 72) of the APRIORI control card, 
internal flags are set (APRKEY^-777777 , APRFLG=APRKEY) to indicate 
that the subclass population data from the statistics file, SAVTAP, 
is to be used for computing subclass apriori probability values. 

If the word FILE is not detected in the parameter field of the 
input APRIORI card, it is assumed that the card contains floating 
point (decimal) numbers, and the parameter field of the card 
is read as is usually done - i.e., scanned for decimal numbers 
separated by a comma, and storing each number found into the 
apriori probability value array. 



Jl tin- 1 laq (APWLG) is sot by virtue of detecting the word FILE 
in the parameter field of the APRIORI control card, subprogram 
SETUP2 obtains the subclass population data from the statistics 
file, SAVTAP via the array KEPPTS . The computation of subclass 
apriori probability values are performed in SETUP2 and the computed 
values are stored in the apriori probability value array, APRIOR. 
The computation of the subclass apriori probability values are 
performed as follows; 

N 



where 

= apriori probability for subclass i 

= number of pixels in subclass i 

K = total number of pixels in all subclasses 

If the word FILE is detected in the parameter field of the 
CATEGORY control card, each class defined on the SAVTAP file 
becomes a category. For each class name there is a correspond int 
category of like name. 

If the word FILE is not detected in the parameter field of the 
CATEGORY card, it is assumed that a category is being defined by 
analyst input, as usual. 

In the subprogram SETUP2 three existing arrays are initialized as 
follows, when CATEGORY FILE card is input; 

KCLSNA = class names from SAVTAP file 

NOCTCL = 1 for number of classes in category I 

CATNAM = class names from SAVTAP file 

and NOCAT = N0CLS2 where NOCAT is the number of categories and 
N0CLS2 is the number of class names on the SAVTAP file. 
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3.2.6 SOFTWARE COMPONENT NO. 6 (DISPLAY) 


The classification results display processor, DISPLAY, is modified 
to meet LACIE Procedure 1 image processing requirements, as 
follows : 

1. to accept a new processor control card 

2. to accept the dot data file, DOTFIL 

3. to provide a dot classification performance summary by dot 
categories which also includes; 

• a tabulation of both the uncorrected proportion and the 
bias corrected proportion of each dot category in the 
total area classified 

• an "alpha" table which tabulates proportions for each 
labeled category of bias correction dots (type 2 dots) , 
which for each category are the ratios of the number of 
dots with the given analyst-labeled category that were 
classified into each possible category, to the total 
number of dots classified into each category. 

4. to provide a dot classification performance summary for each 
dot on the analyst's specified file (DOTFIL). 

3. 2. 6.1 Linkages 

The DISPLAY processor requires the Uni vac EXEC2 operating system, 
the Fortran V compiler, the EOD-LARSYS utility routines, and 
common blocks GLOBAL and DISPL. 

See appendix F for description of modifications to the following 
subprograms: SETUPS, REDIF3, DSPLY2, and DSPLY; for description 

of subprogram TAPLAB; and for description of internal subprograms 
DOTSUM, DTCHK, and DOTPCT. 



3. 2. 6. 2 Interfaces 


As a result of these modifications, the DISPLAY processor 
optionally accepts as input the dot data file, DOTFIL, created 
by the DOTDATA processor. The file is assigned either to Univac 
logical unit F (Fortran unit 8) or to an analyst specified unit. 
(See appendix G for foinnat of file) 

As usual, the DISPLAY processor requires the input of a classifi- 
cation results file, MAPTAP, created by the CLASSIFY processor. 
The file must be assigned to Univac logical unit B (Fortran 
unit 2) (See Appendix D in EOD-LARSYS User Procedures LEC 3984 
for format of file) . 


3. 2. 5. 3 Inputs 

Requires card input described in section 12.4 of User Documentation 
EOD-LARSYS and the following card: 


• NEW CONTROL CARD 


Keyword 

(col. 1) 
DOTFILE 


dot data classification 
performance summaries. 

The parameter, m, designates 
the file to be processed 
by DISOLAY, and n desianates 
the Fortran unit number 
assigned to the input file, 
DOTFTL. 


Parameter (s) 


Function 


(col. 11-72) 
INPUT/UNIT=n,FILE=m 
(Default: UNIT=8, 

FILE=1 

No dot data classi- 
fication performance 
summaries) 


Initiates the input of 
dot data file, DOTFIL, 
from the designated (or 
default, if not designated) 
unit and/or file, and ini- 
tiates the output of the 




• FIELD CARDS 


If the Procedure I option is to be exercised, the only kinds 
of fields that may be input are DO/DU fields (designates 
other or designated unidentifiable) . 

No test fields may be input; the training fields are the bias 
correction dots on the dot data file, DOTFIL. 

The format of the DO/DU field cards and the method of input 
is given in sections 3.1.3 and 12.4.4 of the User Documen- 
tation EOD-LARSYS, LEG 3984. 

3. 2. 6. 4 Outputs 

The DISPLAY processor normally provides an output tape for 
display on the PMIS DAS (Data Analysis Station) , as an option. 

Also, the DISPLAY processor normally provides the following 
output to the line printer: 

• a line-printer map of each classified field on the input 
classification results file (MAPTAP) . 

• classification summaries for each classified field, showing 
the number of pixels classified into and thresholded from 
each subclass, class, and category. 

• classification summaries of training fields, or optionally 
test fields. The summaries are presented by subclass, class, 
and category. 

As a result of the modifications described in section 3.2.6, the 
DISPLAY processor optionally provides dot data classification 
performance summaries (instead of the normal output classification 
summaries described above) if the DOTFILE control card described 
in section 3. 2. 6. 3 is input to the DISPLAY processor. 


The line printer output of dot data performance summaries is as 
follows : 

• A tabulation of individual dots (pixels) showing the sample/ 
line coordinates of each dot, the analyst-labeled category of 
each dot, and the classification category of each dot. 

• A classification performance svuranary by analyst-labeled dot 
data category. This performance summary includes the uncor- 
rected segment (classified field) proportion for each dot 
category, the bias corrected segment (classified field) propor- 
tion for each dot category, and the bias correction "alpha" 
table. In addition, this summary provides the total number of 
dots (pixels) in each analyst-labeled category, the percent 
correct classification for each analyst-labeled category, the 
number of analyst-labeled dots which were thresholded, and the 
number of dots (pixels) of a given analyst-labeled category 
which were classified into the labeled category and classified 
into categories other than the given labeled category. 

3. 2. 6. 5 Storage Requirements 

With the current overlay (MAP) structure, the DISPLAY processor 
requires 16012 locations in the Data Bank, 6732 locations in the 
Instruction Bank, for a total of 32744 locations, plus 19345 
locations for the resident programs. 

3.2.6. 6 Description 

The modifications to the DISPLAY processor are such that the 
analyst may either exercise the normal processor capabilities as 
documented in the User Documentation EOD-LARSYS, LEG 3984, or the 
LACIE Procedure I capabilities. The difference between the two 
capabilities is in the type and format of classification perfor- 
mance tables output. See Section 3. 2. 6. 4 for a description of 
the output tables. 


The processor as modified uses existing available storage, 
normally occupied by training field data, to store the dot 
coordinates and category labels, from the input dot data file. 

The subprograim, REDIF3, is modified to read and de-code the new 
DOTFILE control card, and to set a flag which indicates that 
processing of the dot data is to be performed by DISPIAY. 

Subprogram DSPLY2 is modified to construct the individual dot 
classified category table, using existing storage to construct 
the table, and also to build the dot classification described in 
section 3. 2. 6. 4. The modifications to DSPLY2 are primarily the 
addition of two new internal subprograms. Subroutine DOTPCT, 
which constructs the table of dot classification results, and 
subroutine DOTSUM, which provides the computations and output of 
the two dot classification summaries. 

In subroutine DOTSUM, the computations of the uncorrected category 
proportions, the "alpha" table, and the bias corrected category 
proportions are as follows; 

Pu^ = uncorrected segment (classified field) proportion 
for category 

= number of dots (pixels) classified as 
total number pixels in classified field 


for i = 1 , 2 , . . . ,NOCAT categories of analyst- 

labeled dots. 



Alpha. . “otjj “ Nvunber of dots labeled C., classified as C. 
i/D ij i! 1 

Number of dots classified C^, 

for i,j/ » 1,2,.,., NOCAT categories of labeled 

dots 

Pc^ » bias corrected proportion for category 
AT 

(a. . * Pu- ) 

13 J 

The computation of the percent correct classification of the 
dots is performed as follows: 

Number of dots (pixels) labeled C^, classified as 

Total number labeled dots in 
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3.2.7 SOFTWARE COMPONENT NO. 7 (UTILITY Rcn^TINES) 

Three* EOD-LARSYS system utility subproerrams are describotl in 
this document; TAPERD, RDDOTS, and WRTDOT. 

• TAPERD 

Routine TAPERD reads the MSS data tape (DATAPE) in either the 
Universal or LARSYS 2 format. (Most of the processors in 
the EOD-LARSYS system invoke this routine.) 

• RDDOTS 

Routine RDDOTS reads the dot data tape (DOTFIL) for spectral ' 
spatial information concerning the dots of the Procedure 1 
concept . 

• WRTDOT 

Routine WRTDOT outputs the DOTFIL. In general, two types of 
files are written, the first containing labelina/startina 
dots, and the second containing bias correction dots. 

3. 2. 7.1 Linkages 

The utility subprograms use the Fortran V compiler, Univac soft- 
ware system routines, other EOD-LARSYS utility routines, and 
common blocks GLOBAL and INFORM. See appendix G for description 
of modifications to subprogram TAPERD and for description of 
subprograms RDDOTS and WRTDOT. 

3. 2. 7. 2 Interfaces 

Interfaces with other routines mainly via calling arguments, 
but also may reference common blocks INFORM and GLOBAL. 

3 . 2 . 7 . 3 Inputs 

• TAPERD 

Accepts a multi-file Universal or LARSYS 2 formatted tape. 

The tape may be created either by the EOD-LARSYS system or 
by a computer other than UNIVAC. (See appendix A and B, 
respectively in User Documentation EOD-LARSYS for formats.) 


• RDDOTS 


Accepts a multi- file unformatted Fortran created tape in the 
format defined in appendix H. 

3. 2. 7. 4 Outputs 

• WRTDOT 

Outputs a multi-file unformatted Fortran written tape in the 
format defined in appendix H. 

3. 2. 7. 5 Storage Requirements 

The subprograms TAPERD, RDDOTS, and WRTDOT overlay one another. 
The largest of the 3 routines is TAPERD which requires 8700 
locations of core. 

3. 2. 7. 6 Description 

• TAPERD 

The TAPHDR section is modified to unpack the sun angles from 
the header record of the MSS data tape, DATAPE. Also the 
pixel start number, stop number, and skip factor for line 
and sample is extracted from the header record. 

• RDDOTS 

The subprogram RDDOTS reads the interface tape DOTFIL for 
3 processors, LABEL, ISOCLS, and DISPLAY. Since ISOCLS and 
DISPLAY need different portions of the DOTFIL information, 
a switch, passed via a calling argument, controls which 
portions of the file \ return to the calling routine. 

• WRTDOT 

The subprogram WRTDOT outputs the dot data file, DOTFIL, 
for the processors DOTDATA and LABEL. 



4 . OPERATION 


Batch oriented system, EOD-LARSYS, is operational on the Univac 
1108 Exec 2 located in Building JSC12. Guidelines for job sub- 
mittal are outlined in the IDSD Procedures Manual. Guidelines 
for program procedures for executing Procedure 1 are outlined 
in two documents, section 3 of this document and the User 
Documentation EOD-LARSYS . 




5. TEST PROCEDURE 


To verify the accuracy of the added modifications to the EOD- 
LARSYS system, a segment in Rice, Kansas was processed. 

The entire segment was clustered with the clustering algorithm being 
self-initialized. The resultant unlabeled statistics were labeled 
by the all-of-a-kind procedure. The labeled statistics were then 
used by the sum-of-density classifier to classify the LACIE 
segment. Using bias correction dots (subset of the 209 dots) , 
the classification results were adjusted to allow for analyst 
interpretation of a given dot. 


Note, the shown results are to demonstrate programming accuracy 
only. 

See appendix I for vertif ication run. 
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SUBPROGRAM DESCRIPTIONS 
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TABLE OF CONTENTS 


A.l Software for Subprogram No. l(FLDTYP) 
A. 2 Software for Subprogram No. 2(D0TDAT) 
A. 3 Software for Subprogram No. 3 (DOTS) 

A. 4 Software for Subprogram No. 4(SET13) 


A.l SOF'r^VARE SUBPROGRAM NO. 1 (FLDTYP) 


The subprogram FLDTYP initiates the readina of the fiel<i 
data set defined in section 3. 2. 1.3 and signals the calling 
routine DOTS when all the fields for a given type have been 
processed. 

A. 1.1 LINKAGES 

FLDTYP is called by the routine DOTS and calls the routine 
LAREAD. 

A.l. 2 INTERFACES 

Interfaces with other routines through the common blocks 
INFORM and DOTVEC. 

A.l. 3 INPUTS 

Calling sequence: CALL FLDTYP (FIELDS, STAMNT, $, $, $, IPT, 



VERTEX) 



Parameter 

Diraension 

In/Out 

Definition 

FIELDS 

(4,1) 

Out 

Array for stor- 
ing the field 
data: 

STAMNT 

1 

In 

Parameter for a 
go to Fortran 
statement 

$ 

1 

In 

Exit route if a 
field card is 
encountered 

$ 

1 

In 

Exit route if a 
TYPE card is 
encountered 



Pareuneter 

Dimension 

In/Out 

Definition 

$ 

1 

In 

Exit route if 
a $END* card 
is encountered 

IPT 

1 

In 

Pointer for 
FIELDS array 

VERTEX 

1 

Out 

Array for 
storing the 
field vertices 


A.. 1.4 OUTPUTS 

Stores the type and category names in the common block DOTVEC. 

A. 1.5 STORAGE REQUIREMENTS 
Storage used: Code 212g Data 123g. 

A. 1.6 DESCRIPTION 

The subprogram PLDTYP signals the calling routine when a 
TYPE card, field card or $END* card is encountered by exiting 
through calling arguments. If a TYPE card is encountered, 
the subprogram exits by calling argument number 4; if a field 
card is encountered, the subprogram exits by calling argu- 
ment number 3; if a $END* card is encountered, the subprogram 
exits by calling argument number 5. 

A. 1.7 PLOW CHART 

See figure A-1 for flow chart of the subprogram PLDTYP. 

A. 1.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 





A. 2 SOFTWARE SUBPROGRAM hO. 2 (DOTDAT) 


The subprogr 2 un DOTDAT is the driver routine for the DOTDATA 
processor. 

A. 2.1 LINKAGES 

DOTDAT is called by the loonitor routine MONTOR and calls the 
routines SET13 and DOTS. 

A. 2 . 2 It;:TERFACES 
N/A 

A. 2.3 INPUTS 

Calling sequences CALL DOTDAT (ARRAY , TOP ) 

Parameter Dimension In/Out 

ARRAY 10,600 In 


TOP 1 In 

A. 2 . 4 OUTPUTS 
N/A 


Definition 

Array in blank 
common dimen- 
sioned by the 
parameter top 

TOP = 10,600 


A. 2. 5 STORAGE REQUIREMENTS 
Storage Used: Code - 6g Data - 35g 

A. 2. 6 DESCRIPTION 

The subprogram DOTDAT calls SET13 to read and analyze the 
control cards. DOTS is called to coordinate the functions 
required to output the dot data file, DOTFIL. 



A. 2. 7 PLOW CHART 


See figure A-2 for flow chart of the subprogram DOTDAT. 


A. 2. 8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




DOTDAT 


Figure A-: 



Flow chart for the subprogram DOTUAT 




A. 3 SOFTWARE SUBPROGRAM NO. 3 (DOTS) 

The subprogram DOTS is the coordinator and suoervisor of the 
functions reauired to cutout the dot data file. DOTFIL. This 
file is an interface for the orocessors ISOCLS. LABEL and 
DISPLAY. 

A. 3.1 LINKAGES 

DOTS is called by the routine DOTDAT and calls the routines 
TAPHDR, FLDINT, LINERD, FDLINT, WRTFLD, WRTDOT, and internal 
subprogram PRTDOT. 

A. 3. 2 INTERFACES 


Interfaces with 
GLOBAL, DOTVEC, 

otiier routines through the 
and ISOLNK. 

common blocks INFORM, 

A . 3 . 3 . INPUTS 

Calling sequence 

; CALL DOTS 

(DATA, FIELDS, 

VERTEX, TOP) 

Parameter 

Dimension 

In/Out 

Definition 

DATA 

5000 

In 

Array for 
storing the 
dot data. 

FIELDS 

1000 

In 

Array for stor- 
ing the field 
data 

VERTEX 

1000 

In 

Array for storing 
the field vertices 

TOP 

1 

In 

Parameter. 


TOP=10,600. 




A. 3. 4 OUTPUTS 


Optionally, prints the spatial and spectral information of the 
dot data file, DOTFIL, on the line printer. 

A. 3. 5. STORAGE REQUIREMENTS 

Storage used: Code: 754g Data: 23633g 

,\.3.6 DESCRIPTION 

DOTS calls the routine FLDTYP to coordinate the reading of the 
field cards for a given type. For each field in type i informa- 
tion such as line number, sample number, type number, category 
number, and dot (grid point) is collected. A maximum of 
n£5000/ (number of channels) £250 dots may be collected for type i. 

When all the fields for type i have been processed, the routine 
DOTFIL will output the data to a file. A file is written for 
each type i. 

If applicable, the sun angles from the input imagery data tape 
will be output on the DOTFIL tape. 

A. 3. 7 FLOWCHART 

See figure A-3a for flow chart of the subprogram DOTS and fig- 
ure A-3b for flow chart of the internal subprogram PRTDOT. 

A. 3. 8 LISTING 

Available in the Data Technique Laboratory, JSC Building 17. 


(^1 




i 


0 1 

y 


A 


CALL FDLINT . 
(Find inter- \ 
t.cepts and NI) / 




f- - 

"I 

j J - 0 
J.T - j+2 

1 i 

J > NI i 

I 


J — ■! No 


IB “ be 9 inning 
sample 


I IE “ ending I 

I s^unplc J 


K - IB 


> K "• K+1 


K > IE 

i 


r~ 

TOTVEC “ TOTWC 
+ 1 



Store spatial 
and spectral 
information 
into array 
DATA 


U- 


ORIGINAL PAGE IS 
OF POOH QUALHY 


Figure A-3a.— Flow chart for subprogram DO'rS 
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Figure A- 3b.— Flow chart for internal subprogram PRTDOT 





A. 4 SOFTWARE SUBPROGRAM NO. 4 (SET13) 


The subprogram SET13 reads and analyzes the control cards 
for the DOTDATA processor. 

A. 4.1 LINKAGES 

SET13 is called by DOTDAT and calls the routines NUMBER, 

ORDER, NXTCHR, and FIND. 

A . 4 . 2 INTERFACES 

Interfaces with other routine through the common blocks 
GLOBAL, INFORM, and DOTVEC. 

A. 4. 3 INPUTS 

• Calling sequence: CALL SUBROUTINE SET13 

(No calling arguments are used) 

• Control Cards 

For description of control cards acceptable by this 
routine, see section 3.2.1. 3. 

A . 4 . 4 OUTPUTS 

A line printer summary of the selected options is printed. 

Parameter values such as file switches and option switches 
will be stored in the common blocks. 

A. 4. 5 STORAGE REQUIREMENTS 

Storage used: Code - 62 , Data - 241g. 

A. 4. 6 DESCRIPTION 

Each control card will be read and analyzed for format errors. 
If the control card is valid, the appropriate parameters and 



7 / 


switches will be set. When appliable, if a control card 
for a specific option is not input, a default value will be 
supplied. 

A. 4. 7 FLOWCHARTS 

See figure A- 4 for flow chart of the subprogram SET13. 

A. 4. 8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.l SOFTWARE SUBPROGRAM NO. 1 (ALLKIN) 


The subprogr^un ALLKIN labels the statistics obtained from the 
clustering processor ISOCLS by all-of-a-kind procedure. 

B.1.1 LINKAGES 

The subprogram ALLKIN is called by the routine LABLR and if 
applicable, calls the routine KNEAR. 

B.l. 2 INTERFACES 

Interfaces with other routines through the common blocks INFORM 
and LABS. 

B.l. 3 INPUTS 

Calling sequence; CALL ALLKIN (DOTS, SUBVEC, SUBNO, CATVEC, MEAN) 


Parameter 

Dimension 

In/Out 

Function 

DOTS 

(SIZE, TOTDT2) 

SIZE = 4 + NOFEAT 
TOTDT2 = total dots 
on DOTFIL 

IN 

Array containing 
the dot data 
information 

SUBVEC 

NOSUB2 

NOSUB2=No. of 
clusters 

IN 

Array containing 
the sequence no. 
used for ordering 
the output statis- 
tics 

SUBNO 

NOSUB 2 

OUT 

Array containing 
the no. of clusters 
in each category 

CATVEC 

NOSUB2 

OUT 

Array containing 
the category no. 
each cluster was 
assigned during this 
labeling procedure 


Parameter 

Dimension 

In/Out 

Function 

MEANS 

(N0PET2,1) 

IN 

Array containing 


NOFET2-NO. Of 


the means of the 


channels 


unlabeled statis- 




tics 



B.1.4 OUTPUTS 


A line printer summary of the labeling dots within a cluster 
will be printed. 

B.1.5 STORAGE REQUIREMENTS 

Storage used; Code - 772g Data - 2262g 

B . 1 . 6 DESCRIPTION 

For a given cluster, the labeling dots will be scanned and all 
the dots assigned to that cluster during the clustering process 
will be grouped. 

The labels for the group of dots will be polled. If the labels 
are of one category, the cluster will be labeled that category. 

If more than one category is encountered, the category with the 
majority of the dots will label the cluster. If a given cluster 
does not contain any labeling dots within the cluster, the k- 
nearest neighbor procedure will be applied. 

B . 1 . 7 FLOWCHARTS 

See figure B-1 for flow chart of the subprogram ALLKIN and figure 
B-la for internal subprogram CHECK. 

B.1.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 










Figure li-l.— Flow chart of iubi'rogr<.\n AI.IXIN (continued) 




CHECK 






Discard the 
furtherost 
dot J 


CATGRY(J) = 
CATGRY(J) -1 


Hove dot 
sequence no. 


. 1 ,. 


Find largest 
MAX 



UKi-JIIvAL PAGE) IS 
OF P(’K)T- ? l-T-v 


Figure B-la.— Flow chart for internal s 


iuhluogram CHUCK 



B.2 SOFTWARE SUBPROGRAM NO. 2 (CLS''4AP) 


The subprogram CI.SMAP outputs conditional and/or mixed cluster 
maps. 


B . 2 . 1 LINKAGES 

The subprogram CLSMAP is called by the routine LABLR and calls 
the routines WRTHED, WRTLN, MAPHND, and CLRKEY. 

B . 2 . 2 INTERFACES 

Interfaces with other routines through the common blocks INFORM, 
GLOBAL and LABS. 

B.2. 3 INPUTS 

Calling sequence: C7iLL CLSMAP (CATSUB, SWTCH, SUBNO, SUBVEC , 

SUBDES, CATVEC) 


Parameter 

Dimension 

In/Out 

Definition 

CATSUB 

NOSUB2 

N0SUB2=No. of 
clusters 

In 

Array that flags 
the mixed and 
conditional 
clusters 

SWTCH 

1 

In 

oWTCH=l - condi- 
tional map is 
output 

SWTCH=2 - mixed 
map is output 

SUBNO 

NOCAT 

In 

Array containing 
the no. of clusters 
in a category 

SUBVEC 

NOSUB2 

In 

Array containing 
the sequence no. 
used for ordering 
the output statistics 


SL'BDES 


N0SUB2 


CATVEC N0SUB2 


In Array containing 

the cluster names 

In Array containing 

the category no. 
each cluster was 
assigned 


B . 2 . 4 OUTPUTS 

A conditional and/or cluster map will be output on the analyst 
specified unit. The default unit will be FORTRAN unit 16 (logical 
unit N) . (See section 5.1 in LARSYS document for format of file) . 
If output of both a conditional and mixed cluster map is requested, 
the conditional map is ouhput first; then, the mixed map. 

A line printer symbol map of the conditional and/or mixed cluster 
map will be printed. 

B . 2 . 5 STORAGE REQUIREMENTS 
Code-1034g Data 7360g 

B.2.6 DESCRIPTION 

For the line printer conditional or mixed map, each unconditional 
labeled cluster is assigned a category symbol and each conditional 
or mixed cluster is assigned a unique symbol. 

For the output PCF tape the unconditional labeled clusters are 
arranged in alphabetical order, followed by the conditional or 
mixed clusters, arranged by cluster number in ascending order. 

By a table look-up technique each cluster is assigned a color. 

The previously stored unconditional map (MAPFIL) is read into 
core from drum a line at a time. Each pixel of the line is 




assigned the appropriate predetermined symbol or color and output 
to the line printer and PCF tape, respectively. 

B . 2 . 7 FLOWCHART 

See Figure B-2a and B-2b for flow chart of CLSMAP. 


B.2.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.3 SOFTWARE SUBPROGRT^ NO. 3 (DOTDST) 


The sxibprogreun DOTDST computes the inter-cluster dot distance. 
The distances for all the labeling dots are computed and saved 
on drum for later referencing. 

B.3.1 LINKAGES 

The subprogreun DOTDST is called by the routine LABLR and calls 
the routine SUNFCS. 

B.3.2 INTERFACES 

Interfaces with other routines through the common blocks GLOBAL, 
INFORM, and LABS. 

B.3. 3 INPUTS 

Calling sequence: CALL DOTDST (MEANS, DOTS, TABLE, TOP) 


Parameter 

Dimension 

In/Out 

Definition 

MEANS 

(NOFEAT, NOSUB 2) 

In 

Array containing 
the unlabeled 
means . 

DOTS 

(SIZE, 250) 

In 

Array containing 
the dot data. 

TABLE 

(TOTDT2,l) 

In 

Array containing 
the distances of 
all the dots for 
each cluster. 

TOP 

1 

In 

Parameter that 


defines the maxi- 
miim no. of words 
in TABLE. 





B.3.4 OUTPUTS 


The inter-distances table will be printed. 

B . 3 . 5 STORAGE REQUIREMENTS 
Code-455g Data-504g 

B.3.6 DESCRIPTION 

The distance (LI or L2) from each dot and a cluster is computed 
using a default sun angle correction of 60® or using an analyst 
specified sun angle correction. 

When the distances for a given cluster are computed, they are 
stored on high speed drum a group at a time. (Group = no. of 
labeling dots) . The number of writes to drums depends upon the 
number of cluster means. 

To print the intercluster-dot distance table, the distances are 
read into core, then, printed on the line printer. 

B . 3 . 7 FLOWCHART 

See Figure B-3 for flowchart of DOTDST. 

B,3.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.4 SOFTWARE SUBPROGRAM NO. 4 (KNEAR) 

The subprogram KNEAR labels Srhe statistics obtained from the 
clustering processor ISOCLS by the k-neax'est neighbor procedure. 

B . 4 . 1 LINFJVGES 

The subprogram KNEAR is called by the routine LABLR and does not 
call any LARSYS routines . 

B.4. 2 INTERFACES 

Interfaces with other routines through the common blocks GLOBAL, 
INFORM, and LABS. 

B.4. 3 INPUTS 

Calling sequence; CALL KNEAR (DOTS, SUBVEC, SUBNO, CATVEC, ITER, 



TABl, SWTCH 

, CATNUM, CLUNUM, 

MEANS ) 

Parameter 

Dimension 

In/Out 

Definition 

DOTS 

(SIZE,TOTDT2) 

In 

Array containing 
the dot data in- 
formation 

SUBVEC 

NOSUB2 

In 

Array containing 
the sequence no. 
used for ordering 
the output statis- 
tics . 

SUBNO 

NOSUB2 

Out 

Array containing 
the no. of clusters 
in each category 

CATVEC 

NOSUB2 

Out 

B-17' 

Array containing 
the category no. 
each cluster was 
assigned during tnis 
labeling procedure 



I 


Parameter 

Dimension 

In/Out 

Definition 

ITER 

1 

In 

Number of itera- 
tions 

TA3‘ 

1 

In 

Starting drum 
address for 
retrieving 
distances 

SWTCH 

1 

In 

SWTCH = 0 

CATNUM 

1 

u 

In 

Category number 
cluster CLUNUM 
was assigned 

CLUNUM 

1 

In 

Number of un- 
labeled cluster 
being labeled 

MEANS 

{N0FET2,!) 

In 

Array containing 
the means 


B . 4 . 4 OUTPUTS 

A line printer summary of the k-nearest labeling dots is 
printed. 


B . 4 . 5 STORAGE REQUIREMENTS 

Storage used: Code-613g Data- 1721g 

B.4.6 DESCRIPTION 

For a given cluster, a poll of the category labels of the k-nearest 
dots to the cluster is computed. Th-i label with the majority of 
the dots labels the cluster. If a tie occurs, k-1 dots used. 

The default value of k is 1. 



I 




c 


B . 4 . 7 FLOWCHART 


See figure B-4 for flowchart of KNEAR siibprogram. 

B.4.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 











B.5 SOFTWARE SUBPROGRAM NO. 5 (LABDOT) 


The subprogram LABDOT labels or relabels the aot data file^ 
DOTFIL. 

B.5.1 LINKAGES 

The subprogram LABDOT is called by the routine LABLR and does not 
call any LARSYS routines. 

B . 5 . 2 INTERFACES 

Interfaces with other routines through the common block LABS. 

B.5. 3 INPUTS 

Calling sequence: CALL LABDOT (DOTS) 


Parameter 

Dimension 

In/Out 

Definition 

DOTS 

(SIZE,1) 

In 

Array containing 


the dot data 
information 


B.5. 4 OUTPUTS 
N/A 


B.5. 5 STORAGE REQUIREMENTS 

Storage used: Code-113g Data-120g 


B.5. 6 DESCRIPTION 

The category names input through control cards are checked for 
new name entries. If a new category has been input, the arrays 
concerning the category information is adjusted. 


Each dot number input through a control card in the setup routine 
is relabeled to the appropriate label. No dots may be deleted 




from the set of dot* being updated. 


B . 5 . 7 FLOWCHARTS 

See figure B-5 for flowchart of LABDOT subroutine. 


B.5.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.6 SOFTWARE SUBPROGRAM NO. 6 (LABEL) 


The sxabprogram LABEL is the driver routine for the LABEL processor. 
B.6.1 LINKAGES 

The subprogram LABEL is called by the monitoring routine MONTOR 
and calls routines SET14, LABLR, and PILERD. 

B . 6 . 2 INTERFACES 

Interfaces with other routine through the comiiion blocks INFORM 
and LABS . 

B.6. 3 INPUTS 

Calling sequence: CALL LABEL (ARRAY, TOP) 


Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

10,600 

In 

Array in blank 




common . 

TOP 

1 

In 

Parameter. 




TOP=10,600 


B.6. 4 OUTPUTS 
N/A 

B . 6 . 5 STORAGE REQUIREMENTS 

Storage used: Code-3731g Data-53g 

B.6. 6 DESCRIPTION 

LABEL calls the subprogram SET14 to read and analyze the control 
cards, calls FILERD to read in the specified files, and calls 
LABLR to coordinate the functions required to execute the options 
specified by the analyst. 



B.6.7 FLOWCHART 


See Figure B-6 for flowchart of LABEL. 


B.6.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




Figure B-6.— Flowchart for Subprogram LABEL 



B.7 SOFTWARE SUBPROGRAM NO. 7 (LABLR) 

The subprogram LABLR is the coordinator and supervisor of the 
routines required to perform the operations specified by the 
analyst. 

The analyst may either update the DOTFIL and ■'or SAVTAP files or 
label a given set of statistics by one of two procedures (k- 
nearest neighbor or all-of-a-kind) in one execution of the LABEL 
processor. 

B.7.1 LINKAGES 

The subprogram LABLR is called by the routine LABEL and calls 
routines LABMAN, LABDOT, KNEAR, ALLKIN, CNDMAP, MIXMAP, MANORD, 
REDDER, WRTDOT, DOTDST, MSTA, and DSPTAP. 

B.7. 2 INTERFACES 

Interfaces with routines through the common blocks INFORM, GLOBAL, 
and LABS. 

B.7. 3 INPUTS 

Calling sequence; CALL LABLR (ARRAY, TOP, NOFLD, TOTVRT, FLDSAV, 



VERTX, MEANS) 



Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

10,600 

In 

Array in blank 
common 

TOP 

1 

In 

Parameter. 

TOP=10,600 

NOFLD 

1 

In 

No. of fields on 
the DOTFIL 

TOTVRT 

1 

In 

No. of vertices 
of fields on DOTFIL 



VERTX 


(2,T0TVRT) 


MEANS (N0FET2,1) 


In Array containing 

the field vertices 
of fields on 
DOTFIL 

In Array for storing 

the means 


B . 7 . 4 OUTPUTS 


N/A 


B . 7 . 5 STORAGE REQUIREMENTS 

Storage used: Code-1117g Data-7504g 

B.7.6 description 

Depending on the analyst selected options LABLR supervises the 
following functions: 

• Relabeling of the SAVTAP file 

• Relabeling of the DOTFIL file 

• Execution of the k-nearest neighbor procedure 

• Execution of the all-of-a-kind procedure 

• Output of the conditional cluster map 

• Output of the mixed cluster map 

• Output of the DISPLAY interface tape, MAPTAP 

B . 7 . 7 FLOWCHART 

See figure B-7a and B-7b for flowchart of LABLR subprogram. 
B.7.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.8 software subprogram no. 8 (MANORD) 

The subprogram LABMAN relabels the statistics file, SAVTAP^ by 
input of the control card STATLABEL. 

B.8.1 LINKAGES 

The subprogram MANORD is called by the routine LABLR and does not 
call any LARSYS routines. 

B.8.2 interfaces 

Interfaces with other routines through the common block LABS. 

B.8. 3 INPUTS 

Calling sequence: CALL LABMAN (CLSNAM, CLSVEC, SUBVEC, N0CLS2, 

SUBNO, N0SUB2) 


Parameter 

Dimension 

In/Out 

Dimension 

CLSNAM 

NOCLS2 

In 

Array containing 
class names 
(category names) 
on SAVTAP file. 

CLSVEC 

NOSUB2 

In 

Array containing 
the class-subclass 
correspondence . 

SUBVEC 

N0SUB2 

Out 

Array containing 
the sequence no. 
used for ordering 
the output statis- 
tics . 

N0CLS2 

1 

In 

No. of class 
names on SAVTAP 
file . 




SUBNO 


60 


Out 


Array containing 
the no. of clusters 
in each class 
(category) . 

NOSUB2 1 In No. of clusters 

B . 8 . 4 OUTPUTS 
N/A 

B . 8 . 5 STORAGE REQUIREMENTS 
Storage used: Code-235g Data-166g 

B.8.6 DESCRIPTION 

The analyst may manually relabel the statistics file by input 
of control cards. The classnames (categories) input must match 
a name on the statistics file. 

The analyst may regroup any clusters into another class. For 
example, if the cluster WHTl belonged to the class WHEAT and the 
analyst decided that the cluster should belong to the NONWHT class, 
he may force the cluster WHTl to be assigned to the class NONWHT. 

As a result of regrouping the clusters, the order of the statistics 
on the file are rearranged, but the regrouped cluster are 
not be renamed. 

B . 8 . 7 FLOWCHART 

See figure B-8 for flowchart of MANORD subprogram. 

B.8.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 





















B.9 


SOFTWARE SUBPROGRAM NO. 9 (MIXMAP) 


The subprogreun MIXMAP flags the conditional clusters. 

B . 9 . 1 LINKAGES 

The subprogram MIXMAP is called by the routine LABLR and does not 
call any LARSYS routines. 

B . 9 . 2 INTERFACES 

Interfaces with other routines through the common blocks GLOBAL 
and LABS. 


B.9. 3 INPUTS 


Calling sequence; CALL MIXMAP (DOTS, MIXSUB, NOSUB2 , CATVEC) 


Parameter 


Dimension 


In/Out 


DOTS 


(SIZE,T0TDT2) 


In 


MIXSUB 


N0SUB2 


Out 


NOSUB2 

CATVEC 


1 

NOSUB2 


In 

In 


B.9. 4 OUTPUTS 
N/A 

B . 9 . 5 STORAGE REQUIREMENTS 
Code-156g Data-427g 


Definition 

Array containing 
the dots. 

Array used for 
flagging a mixed 
cluster on the 
output mixed map. 

No. of clusters. 

Array containing 
the category no. 
of each cluster. 



B.9.6 DESCRIPTION 


The cluster number assigned to each dot during clustering is 
retrieved from high speed drum. All the dots belonging to 
cluster I are collected, and their labels are polled. If the 
labels define more than one category, cluster I is flagged as a 
mixed clus _er. A mixed cluster is assigned a cluster of 62-J 
where J = no. of mixed clusters encountered at that time. 

B . 9 . 7 FLOWCHART 

See Figure B-9 for flowchart of MIXMAP. 

B.9.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 











B.IO SOFTWARE SUBPRCXSRAM NO. 10 (SET14) 


The subprogram SET14 reads and analyzes the control cards for 
the LABEL processor. 

B.10.1 LINKAGES 

The subprogram SET14 is called by the driver routine LABEL 
and calls the routines NUMBER, ORDER, NXTCHR, FIND, CRDSCN, 
FLTNUM, and CRDSTA. 

B.IO. 2 INTERFACES 

Interfaces with other routines through the common blocks 
GLOBAL, INFORM, and LABS. 

B.IO. 3 INPUTS 

• Calling sequence: CALL SET14 (ARRAY,TOP) 


Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

10,600 

In 

Array in blank 




common. 

TOP 

1 

In 

Parameter . TOP=l 0,600 


• Control Cards 

For description of control cards acceptable by this routine, 
see section 3. 2. 2. 3. 

B.IO. 4 OUTPUTS 

A line printer summary of the selected options is printed. 

Parameter values such as file switches and option switches 
are stored in the common blocks INFORM, GLOBAL, and 
LABS . 




B . 1 0 . 5 STORAGE REQU I REMENTS 


Storage used: Code-2550g Data-575g 

B.10.6 DESCRIPTION 

Each control card is read and analyzed for format errors. 

If the control card is valid, the appropriate parameters and 
switches are set. When applicable, if a control card for 
a specific option is not input, a default value is supplied. 

B.10.7 FLOWCHART 

See figure B-10 for flowchart of SET14 subprogr 2 im. 


B.10.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




B.ll SOFTWARE SUBPROGRAM NO. 11 (MAPHND) 


The subprogram MAPHND prints out the heading for the conditional 
and mixed cluster maps. 

B.xl.l LINKAGES 

The subprogram MAPHND is called by the routine CLSMAP. 


B.ll. 2 INTERFACES 

Interfaces with other routines through the common block INFORM. 


B.ll. 3 INPUTS 


Calling sequence; CALL MAPHND (NOCAT, CLSSYM, CATNAM, KATNO, 

SUBDES, CATSUB) 


Parameter 

NOCAT 

CLSSYM 

CATNAM 


KATNO 


SUBDES 


Dimension 

1 

NOSUB2 

NOCAT 


N0SUB2 


N0SUB2 


In/Out 

Function 

In 

No. of categories 

In 

Array containing 
the symbols . 

In 

Array containing 
the category 
names . 

Array containing 
the category no. 
each cluster was 
assigned . 

In 

Array containing 
the cluster names 

In 

Array containing 
the flagged 
conditional or 
mixed clusters 


//(■ 


CATSUB 


NOSUB2 


B.11.4 OUTPUTS 


N/A 


B.11.5 STORAGE REQUIREMENTS 
Storage used; Code-175g Data-130g 

B.11.6 DESCRIPTION 

All the clusters for a category are grouped together. Information 
such as the category neune, cluster number, cluster name, and 
symbol for the cluster is printed on the heading of the condi- 
tional and/or mixed cluster map. 

B.11.7 FLOWCHART 

See Figure B-11 for flowchart of MAPHNH iubprogram. 

B.11.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




B.12 SOFTWARE SUBPROGRAM NO. 12 (ASCEND) 


The subprogreun ASCEND sorts an array of floating point numbers 
in descending order. 

B.12.1 LINKAGES 

The subprogreun ASCEND is called by the routine KNEAR and does 
not call any LARSYS routines. 

D.12.2 INTERFACES 

-erfaces with other routines via calling arguements. 

B.12. 3 INPUTS 


Calling sequence: 

CALL ASCEND 

(SCN, LNCAT, 

PTRl, PTR2) 

Parameter 

Dimension 

In/Out 

Def inxtion 

SCN 

LNCAT 

In/Out 

Array to be 
sorted. 

LNCAT 

1 

In 

No. of elements 




to sort. 

PTRl 

LNCAT 

In/Out 

Array containing 
information 
relating to SCN. 

PTR2 

LNCAT 

In/Out 

Array containing 


information 
relating to SCN. 


B.12.4 OUTPUTS 


N/A 


B.12. 5 STORAGE REQUIREMENTS 
Storage used: Code-174g Data-27g 


- ' page io 

i Oi'Aury 


B.12.6 DESCRIPTION 

The subprogram ASCEND operates on 3 arrays, a floating point array 
and 2 integer arrays. The floating point array is arranged in 
decreasing order; the 2 integer arrays are sorted based on the 
ordering of the floating point array. 

B.12.7 FLOWCHART 

See Figure B-7 for flowchart of ASCEND subprogram. 

B.12.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.13 SOFTWARE SUBPROGRAM NO. 13 (CRDSCN) 

The function CRDSCN interprets the control cards DOTLABEL and 
STATLABEL . 

B.13.1 LINKAGES 

The subprogram CRDSCN is called by the routine SET14 and calls 
the routines NXTCHR, FIND, and NUMBER. 

B.13. 2 INTERFACES 
N/A 

B.13. 3 INPUTS 

Calling sequence; CALL CRDSCN (CARD, GRPDEX, GRPNAM, GROUPS, 

NOGRP, GRPTR) 


Parameter 

Dimension 

In/Out 

Definition 

CARD 

62 

In 

Array containing 
the string of 
characters from 
control card 
DOTLABEL or 
STATLABEL . 

GRPDEX 

60 

Out 

Array containing 
the pointers for 
retrieving infor- 
mation from GROUPS 
array. 

GRPNAM 

60 

Out 

Array containing 
the input cate- 
gory names . 




GROUPS 


120 


Out 


Array containing 
the no. of clusters 
associated with 
an input category 
name, followed 
by the actual 
cluster no. 

NOGRP 60 Out No. of category 

names encountered 
on input card at 
a given time. 

GRPTR 1 In Pointer for storing 

into the GROUPS 
array . 


B.13.4 OUTPUTS 
N/A 


B.13.5 STORAGE REQUIRED 

Storage used: Code-345g Data~732g 


B.13.6 DESCRIPTION 

The function CRDSCN interprets an alphanumeric string, one 
character per word. The first 1 to 7 characters are assumed to 
be alpha characters, followed by at least one niomeric character. 

The alpha characters define a category name plus a comma, the 
numeric characters, cluster number (s). 

Three arrays are output. One array contains the labels (category 
names) ; one, the no. of the clusters or dot to relabel; and the 
other, index pointers. 


For the following control cards 




DOTLABEL WHEAT, 3, 4, 7 

DOTLABEL NONWHT, 10, 12 

the array will contain: 

GRPNAM(l) = WHEAT GRPDEX(l). = 1 GROUPS (1) =3 

GRPNAM(2) = NONWHT GPPDEX(2) = 5 GROUPS (2) =3 

GROUPS (3) =4 
GROUPS (4) =7 
GROUPS (5) =2 
GROUPS (6) = 10 
GROUPS (7) = 12 

B.13.7 FLOWCHART 

See figure B-13 for flowchart of CRDSCN subprogram. 

B.13.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.14 SOFTWARE SUBPROGRAM NO. 14 (FILERD) 


The subprogram FILERD computes the needed high speed drum addresses 
and reads in the input files. 

B.14.1 LINKAGES 

The subprogram FILERD is called by the driver routine LABEL and 
calls the routines LAREAD, STOMAP, REDSAV, RDDOTS, and WRTFLD. 

B.14. 2 INTERFACES 

Interfaces with other routines through the common blocks INFORM, 
GLOBAL, and LABS. 


B.14. 3 INPUTS 


Calling sequence: CALL FILERD 

VERTX) 

Parameter Dimension 

ARRAY 10600 

area and contains 
the dot data 
information . 


(ARRAY, TOP, NOFLD, TOTVRT, FLDSAV, 

In/Out Definition 

Out Used for scratch 


TOP 

NOFLD 

TOTVRT 

FLDSAV 


1 

In 

Parameter TOP=10 

1 

Out 

No. of fields on 
DOTFIL. 

1 

Out 

No. of field 
vertices on DOTFIL 

(4, NOFLD) 

Out 

Array containing 
the field informa- 
tion on DOTFIL. 


Out Array containing 

the field vertices 
on DOTFIL. 


2^4 


VERTX 


(2, TOTVRT) 


B.14.4 OUTPUTS 


N/A 

B . 14 . 5 STORAGE REQUIREMENTS 

Storage used; Code-304g Data-105g 

B.14.6 DESCRIPTION 

If a MAPFIL file is input, the subprogram FILERD initiates the 
reading of the field card, computes the dimensions of the input 
file, reads in the file, and stores the entire image on high 
speed drum. 

If a SAVTAP file is input, the statistics are read into the upper 
partition of the array. Only the means and covariances are 
stored on high speed drum. 

If a DOTFIL file is input, the file is read in and stored starting 
in the location of ARRAY previously occupied by the covariances. 

Based on the number of files read in, drum addresses that will 
be needed in other routes are computed and stored in the common 
blocJ: LABS . 

B.14.7 FLOWCHART 

See figure B-14 for flowchart of FILERD subprogram. 

B.14.8 LISTING 

Available in the Data Technique Laboratory, JSC Building 17. 









B.15 software subprogram no. 15 (STOMAP) 


The subprogram STOMAP reads the input MAPFIL tape and stores the 
image on high speed drum. 

B.15.1 LINKAGES 

The subprogram STOMAP is called by the routine FILERD and calls 
the routines TAPHDR, FLDINT and LINERD. 

B.15. 2 INTERFACES 

Interfaces with other routines through the common block GLOBAL. 

B.15, 3 INPUTS 

• MAPFIL tape 

Accepts the cluster or classified map output by ISOCLS or 
DISPLAY. See section 5.1 of User Documentation EOD-LARSYS. 

• Calling arguments: CALL STOMAP (ILINE, NSAMP, HIST, LIMIT, 

BEGINl) 


Parameter 

Dimension 

In/Out 

Definition 

ILINE 

1 

In 

No. of lines on 
MAPFIL 

NSAMP 

1 

In 

No. of samples 
per scan line. 

HIST 

LIMIT 

In 

Scratch area 

LIMIT 

1 

In 

Parameter 
dimensioning HIST 

BEGINl 

1 

In 

Beginning drum 
address for stor- 
ing MAPFIL 




B.15.4 OUTPUTS 
N/A 


B.15.5 STORAGE REQUIREMENTS 
Storage used; Code-241g Data-57g 

B.15.6 DESCRIPTION 

Based on the size of both the array HIST, and the input file the 
subprogram STOMAP computes how many scan lines of the MAPFIL file 
can be read into core at one time. When the buffer is filled, 
the data is dumped to high speed drum. (The number of dumps to 
the drum is a function of the number of lines that can reside in 
core) . 


B.15.7 FLOWCHART 

See Figure B-15 for flowchart of subprogram STOMAP. 

B.15.8 Listing 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B.16 SOFTWARE SUBPROGRAM NO. 16 (DSPTAP) 


Subprogram DSPTAP writes an unformatted unconditional cluster map 
tape to be used as input to the DISPLAY processor to obtain certain 
general performance summaries. The layout and structure of this 
tape is described in the User Documentation EOD-LARSYS manual 
in the SYSTEM INTERNAL FILES section (MAPTAP) . It has 
the same layout as the classification map tape normally 
written by the CLASSIFY processor. However, only 
category/subclass information is written. 

The trigger to initiate writing of this tape is set by inclusion 
of the MAPTAP control card as input to the LABEL processor. 

B.16.1 LINKAGES 

Subprogram DSPTAP is called by routine LABLR. DSPTAP calls the 
UNIVAC tape routine FSBSFL and the drum read routine RREAD. 

B.16. 2 INTERFACES 

Interfaces with other routines through common blocks INFORM, 

GLOBAL, and LABS. 

B.16. 3 INPUTS 


Calling sequence; 

DSPTAP (SUBNO, 
SUBVEC, MEANS, 

SUBDES, FLDSAV, VERTX, CATVEC , 
COVAR, TOP, DATA, NOFLD, TOTVRT) 

Parameter 

Dimension 

In/Out 

Definition 

SUBNO 

NOCAT 

In 

Number of sub- 
classes in each 
category . 

SUBDES 

N0SUB2 

In 

Subclass names 




FLDSAV (4,N0FLD) In For each training 

field 





1. field nzime 

2. class no. of 
field 

3. subclass no, 
of field 

4. no. of field 
vertices 

VERTX 

(2,TOTVRT) 

In 

Field vertices 

CATVEC 

60 

In 

Category no for 
each subclass 

SUBVEC 

60 

In 

Re-ordering of 
original subclass 
order by category 
niomber 

MEANS 

(NOFET2,NOSUB2) 

In 

Subclass mean 
vectors 

COVAR 

(VARSZ2, NOSUB 2) 

In 

Subclass covariance 
matrices, by row 
as lower triangular 

TOP 

1 

In 

Upper limit of 
ARRAY vector, used 
as general buffer 
storage. TOP set 
in MONTOR routine, 
currently 10600 

DATA 

(PTS, NOLINE) 
NOLINE=No. of 
lines real 
into core 

In 

Original cluster 
number the pixel 
was assigned 
during clusterina 
(in ISOCLS) 



NOFLD 1 In Total no. of 

fields 

TOTVRT 1 In Total no. of 

vertices 

B.16.4 OUTPUT 

An unformatted tape (MAPTAP) is written. 

B . I 'S . 5 STORAGE REQUIREMENTS 
Code-1005g Data-4223g 

B.16.6 DESCRIPTION 

An unformatted unconditional cluster map tape is written for use 
by the DISPLAY processor. The layout is the same as that for the 
MAPTAP written by the CLASSIFY processor. However, only category 
and subclass information is provided, there is no inclusion of 
class/subclass and category/class associations. Thxe information 
is that obtained as a result of cluster labeling in the LABEL 
processor, either by tie k-nearest neighbor or all-of-a-kind 
criteria . 

B.16.7 FLOWCHART 

See Figure B-16 for flowchart of DSPTAP subprogram. 

B.16.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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Figure B-16.— Flow chart for subprogram DSPTAP 






B.17 SOFTWARE SUBPROGRAM NO. 17 (NAMSTA) 


Subprogram NAMSTA provides cluster names based on labeling results 
in the LABEL processor. These names are based on category names 
and seciuence numbering. A check is made to ensure that each 
category is involved in the naming of at least one cluster. If 
not, that category is deleted from the total count of number of 
categories . 

B.17.1 LINKAGES 

Subprogram NAMSTA is called by routine LABLR. 

B.17. 2 INTERFACES 

Interfaces with other routines through common blocks INFORM and 
LABS- 


B.17.1 INPUTS 

Callincj sequence: NAMSTA (SUBNAM, 

Parame ter Dimension 

SUBNAM 60 

SUBVEC 60 


N 60 

CATVEC 6 C 


SUBVEC, N, CATVEC, SUBNO) 


Tn/Out 

Definition 

Out 

Cluster (subclass) 
names 

In 

Not used in 
NAIVISTA, re-order- 
ing of original 
subclass names by 
category 

In 

Not used 

In 

Category no. of 
the subclasses 

In 

No. of subclasses 
in each category 


B-'6l 


SUBNO 


60 


B.17.4 output 


N/A 


B . 17 . 5 STORAGE REQUIREMENTS 
Code-203g Data 40g 

B . 17 . 6 DESCRIPTION 

Cluster names are stored in vector SUBNAM on the basis of the 
category assignment to each cluster. These names are six 
characters long. The first three characters are derived from 
the category name, and the last three correspond to sequence 
numbers, starting with '001' for each category. In the event 
that a category has no associated clusters, SUBNO and NOCAT are 
adjusted to remove that category from subsequent consideration. 

B . 17 . 7 FLOWCHART 

See Figure B-17 for flowchart of NAMSTA. 

B.17.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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Figure B-17.— Flow chart for subprogram NAMSTA 
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B.18 SOFTWARE SUBPROGRAM NO. 18 (REODER) 


Subprogram REODER provides a re-ordering of subclass names and 
number of associated pixels, as part of the manual relabeling 
of stats option in the LABEL processor. 

B.I8.1 LINKAGES 

Subprogram REODER is called by routine LABLR. 

B.18. 2 INTERFACES 

Interfaces with other routines through the common block INFORM. 
B.18. 3 INPUTS 


Calling sequence: 

REODER (ARRAY, 

SUBVEC, N) 


Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

TOP (currently 
10600) 

In/Out 

Subclass names, 
starting at index 
SUBDS2 

SUBVEC 

60 

In 

Reordering of 
subclass by 
category 

N 

60 

Out 

No. of pixels in 
each subclass, 
reordered using 
SUBVEC 


B.18. 4 OUTPUTS 
N/A 


B.18. 5 STORAGE REQUIREMENTS 
Code-107g Data-124g 




B.18.6 DESCRIPTION 


The subclass names and number of pixels per subclass are reordered 
according to SUBVEC. The new ordering places all subclasses 
associated with the first category first in order, all subclasses 
associated with the second category second in order, and so on. 

B.18.7 FLOWCHART 

See Figure B-18 for flowchart of the REODER subprogram. 

B.18.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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Figure B-18.- Flow chart for subprogram Ri;ODER 




B.19 SOFTWARE SUBPROGRAM NO. 19 (LABMAN) 


Subprogram LABMAN, under the manual relabeling of stats option 
of the LABEL processor, writes a statistics tape (SAVTAP) and 
optionally punches a card deck. 

B.19.1 LINKAGES 

Subprogram LABMAN is called by the LABLR routine. LABMAN calls 
the UNIVAC tape file positioning on two occasions, one for manual 
stat relabeling, the other for the LABEL routine FSBSFL and the 
drum read routine RREAD. 

B.19. 2 INTERFACES 

Interfaces with no routines other than LABLR, and that only 
through the calling sequence. 

B.19. 3 INPUTS 

Calling sequence: LABMAN (UNIT, FILE, NOCLS, TOTSUB, NOFEAT, 

TOTFLD, TOTVRT, FETVEC, FLDSAV, VERTEX, CLSNMS , NOSUBS, SUBNM, N, 
STADRS, VARSIZ, PUNCH, SUBVEC, PRNSTS) 


Parameter 

Dimension 

In/Out 

Definition 

UNIT 

1 

In 

Statistics tape 
unit number 

FILE 

1 

In 

Statistics tape 
file no. for this 
output file 

NOCLS 

1 

In 

No. of classes 

TOTSUB 

1 

In 

No. of subclasses 

NOFEAT 

1 

In 

No. of channels 
(features) 

TOTFLD 

1 

In 

No. of fields 






\ 

i 

! 


FETVEC 

NOFEAT 

In 

Features | 

FLDSAV 

NOFLD2*4 

In 

Field information 1 

! 

by field (name, j 

class, dummy | 

entry, no. of 
vertices in 
fields 

CLSNMS 

NOCLS2 

In 

Class names 

NOSUBS 

1 

In 

No. of clusters 
in each class 

SUBNM 

NOSUB 2 

In 

Reordered subclass 
names 

N 

60 

In 

Reordered no. of 
pixels in each 
subclass 

STADRS 

1 

In 

Starting address 
for drum storage 
access 

VARS I Z 

1 

In 

Storage size for 
lower triangular 
part of covariance 
matrix 

PUNCH 

1 

In 

Trigger to punch 
card deck 

SUBVEC 

60 

In 

Reordering of 
subclasses by 
category 

PRNSTS 

1 

In 

Trigger to print 
statistics sumntarv 


/¥i 


B.19.4 OUTPUT 


An unformatted statistics tape (SAVTAP) is output. On 
options, a punch card deck and/or a line printer summary is 
performed . 

B . 19 . 5 STORAGE REQUIREMENTS 
Code-1222g Data-1350g 

B.19.6 DESCRIPTION 

A statistics tape is written under the SAVTAP format to reflect 
the manual relabeling of clusters under this option of the LABEL 
processor. On option, a punched card deck is produced. 

B.19.7 FLOWCHART 

See Figure B-19 for flowchart of LABMAN. 

B.19.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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Figure B-19.— Flow chart for subprogram LABKAN 







B.20 SOFTWARE SUBPROGRAM NO. 20 (SUNFCS) 

Subprogrzun SUNFCS computes sun angle gain corrections for pixel 
radiance values for each channel based on input sun angles. 

B.20.1 LINKAGES 

Subprogram SUNFCS is called by routine DOTDST. 

B.20. 2 INTERFACES 

No common blocks are provided. 

B.20. 3 INPUTS 


Calling sequence: 

SUNFCS (SUNCOR, 
ISUNT) 

SUNANG, 

FETVEC, NOFEAT, ISUNC, 

Parameter 

Dimension 

In/Out 

Definition 

SUNCOR 

30 

Out 

Corrections 
corresponding to 
channels in FETVEC 

SUNANG 

8 

In 

Sun angles 

FETVEC 

30 

In 

Channels 

NOFEAT 

1 

In 

No. of channels 

ISUNC 

1 

In 

Flag indicating 
sun angles came 
from cards. 

I SUNT 

1 

In 

Flag indicating 
sun angles were 
extracted from 
MSS data tape 


B.20. 4 OUTPUT 
N/A 


/ 


B . 2 0 . 5 STORAGE REQU I REMENTS 


Code-254g Data-576g 

B.20.6 DESCRIPTION 

The sun angle corrections are computed from an internal table 
covering four channels at a haze factor of 25 percent. The 
corrections correspond to the channels contained in FETVEC. If 
the ISUNT flag is on, the first sun angle in SUNANG is assumed to 
correspond to pass 1. If the ISUNC flag is on, the first sun 
angle in SUNANG is assumed to correspond to the first pass in 
FETVEC . 


The sun angle corrections, which are basica?.ly gain factors, are 
used in the computation of cluster mean-dot distances. 

B.20.7 FLOWCHART 

See Figure B-20 for flowchart of SUNFCS. 

B.20.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 








B.21 SOFTWARE SUBPROGRAM NO. 21 (CLRKEY) 


Subprogram CLRKEY writes the color keys on the mixed and condi- 
tional cluster map tapes. 

B.21.1 LINKAGES 

Subroutine CLRKEY is called by routine CLSMAP. CLRKEY calls 
WRTLN, an entry point in TAPWRl . 

B.21.2 INTERFACES 

There are no common blocks associated with CLRKEY. 

B.21. 3 INPUTS 


Calling sequence: 

CLRKEY (XSIZ, 

I DATA, N0SUB2, 

CH, MEANS, NC) 

Parameter 

Dimension 

In/Out 

Definition 

XSIZ 

1 

In 

No. of pixels/ 
line (maximum 
200) 

I DATA 

(XSIZ,CH) 

Out 

Colors output 
per line 

NOSUB2 

1 

In 

No. of subclasses 

CH 

1 

In 

Set equal to 1; 
in the multi- 
channel case, set 
equal to NC+1 

MEANS 

(NC,NOSUR2) 

In 

Colors corresponding 
to subclasses. 

NC 

1 

In 

Sot equal to 1, 
the no. of channels 




/vf 


B.21.4 OUTPUT 


The color keys are added to the mixed or conditional map tape . 
(See appendix H for format.) 

B.21.5 STORAGE REQUIREMENTS 
Code-247g Data-51g 

B.21.6 DESCRIPTION 

For each cluster (subclass) , a 10 x 10 color number is written 
on the conditional and/or mixed cluster map tapes. This cluster 
number is obtained from an internal table. 

B.21.7 FLOWCHART 

See Figure B-21 for flowchart of CLRKEY program. 

B.21.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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B . 2 2 S OFTWARE SUBPROGRAM NO. 22 (CNDMAP ) 


Subprogram CNDMAP flags conditional clusters. 

B.22.1 LINKAGES 

Subprogram CNDMAP is called by routine LABLR. CNDMAP calls the 
UNIVAC drum read routine RREAD. 

B.22.2 INTERFACES 

Interfaces with other routines through common blocks INFORM and 
LABS. 


B.22.3 INPUTS 


Calling sequence: 

CNDMAP (DOTS, 

CNDSUB, CATVEC) 


Parameter 

Dimension 

In/Out 

Definition 

DOTS 

( 4 , TOTDT2 ) 

In 

Dot data-for each 
dot 

1. sample number 

2. line number 

3 . type 

4 . category no 

CNDSUB 

60 

Out 

For each subclass: 
set equal to entry 
in CATVEC if not 
conditional, other 
wise set equal to 
62 minus no. of 
conditional sub- 
classes flagged 
previously 

CATVEC 

NOSUB2 

In 

Category no. of 
each subclass 




B.22.4 OUTPUT 


N/A 

B.22.5 storage REQUIREMENTS 
Code-142g Data-424g 

B.22.6 DESCRIPTION 

For each cluster I, the subprogram CNDMAP reads into core from 
drum that portion of the distance table for cluster T. The 
minimum distance between the labeled cluster I and dot of like 
label is found and compared with a threshold value. If the dis- 
tance is greater than the threshold, the cluster is flagged as 
conditional. This flagging is done by assigning it a cateogory 
number equal to 62 minus the number of clusters already flagged. 

B.22.7 FLOWCHART 

See Figure B-22 for flowchart of CNDMAP. 

B.22.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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APPENDIX C 


SUBPROGRAM DESCRIPTIONS 


APPENDIX C 


TABLE OF CONTENTS 
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1 (ISODAT) 

2 (RDDATA) 

3 (SETUP?) 

4 (SUNFAC) 

5 (PSPLIT) 


C.l SOFTWARE SUBPROGRAM HO. 1 (ISODAT) 


The new capeQ^ilities of subprogram ISODAT are: 

1. on option, call the utility routine RDDOTS to obtain 
cluster starting vectors, 

2. on option, apply a svm angle correction by channel to 
each pixel during clustering, 

3. remove "designated other" or "designated unidentifiable" 
pixels from the clustering algorithm input. 

C.1.1 LINKAGES 

The Subprogram ISODAT is called by the routine ISOCLS, and 
calls routines PRINT, PSPLIT, CLDIST, and DELETE. 

C.1.2 INTERFACES 

Subprograun ISODAT interfaces with other routines through 
common blocks GLOBAL, PASS, and JSOLNK. 

C.l. 3 INPUTS 

• Calling sequence: CALL ISODAT (C, IPLACE, MEANS, N, STDEV, 

CLD, FLDINF, AVP, AMN) 


Parameter 

Dimension 

In/Out 

Definition 

C 

( NOFEAT, NOPTS) 

Out 

Storage for one 
drum read of pixels 
for clustering 

IPLACE 

NOPTS 

Out 

NOPTS is the number 
of points for one 
drum write of pixel 
cluster numbers. 



5 


Parameter 

Dimension 

In/Out 

Definition i 

j 

MEANS 

(NOFEAT ,MAXCLS) 

In/Out 

Starting mean values ! 

for clustering, MAXCLS 
is the maximum number 
of clusters. 

N 

MAXCLS 

Out 

Number of pixels per 
cluster. 

STDEV 

( NOFEAT, MAXCLS) 

Out 

Standard deviations 
per cluster/channel. 

OLD 

(I lAXCLS, MAXCLS) 

In 

Cluster distance 
matrix. 

FLDINF 

(IPT) 

out 

Field inforraation 
(vertices, vertex 
location) . 

AVP 

(NOFEAT, MAXCLS) 

Out 

Temporary storage for 
computation of cluster 
means 

AMN 

(NOFEAT, MAXCLS) 

In 

Temporary storage for 
computation of cluster 
standard deviations. 


C . 1 . 4 OUTPUTS 
No additional outputs . 

C.1.5 STORAGE REQUIREMENTS 
Code-2176g Data-4400g 


C;^ 


C.1.6 DESCRIPTION 


A dot file trigger in labeled comraon ISOLNK invokes a call to 
RDDOTS for starting vectors. The unit and file numbers will 
have been loaded in labeled common GLOEAL. Sun angle correction 
triggers in labeled common ISOLNK causes the sun angles (stored 
in ISOLNK) to be used to find sun angle corrections from an 
internal table. Application of the correction consists of 
multiplying the pixel channel radiance values by the corrections 
in the pixel assignment to clusters, and the computation of 
maximum standard deviations. 

The DO/DU pixels are identified from the input DO/DU field 
information. As regular field pixels are extracted from the MSS 
data tape in routine RDDATA, any DO/DU pixels contained therein 
are identified and the radiance valves are set to 0 or 255 for 
subsequent identification and deletion in the clustering procedure. 

C . 1 . 7 FLOWCHART 
N/A 

C.1.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




C.2 SOFTWARE SUBPROGRAM NO. 2 (RDDATA) 


Tne new capabilities of subprogram RDDATA are: 

1. read DESIGNATED OTHER and DESIGNATED UNIDENTIFIABLE class- 
n^une cards and the associated field definition cards, and 
calculate tne field information (number of vertices, vertex 
locations, nximber of pixels) for both classes, 

2- compute information concerning the DO/DU pixels fur use 
by other subprograms. This includes the number of DO/DU 
pixels in each pixel field for the current class, and 
drum addresses for DO/DU pixies, 

3. assign printing symbols to DO and DU pixels, as well as 
special mean values (0 or 255) . 

C.2.1 LINKAGES 

The subprogram RDDATA is called by the routine ISOCLS, and 

calls the routines TAPERD and LAREAD. 

C.2. 2 INTERFACES 

Subprogram RDDATA interfaces with other routines through com- 
mon blocks GLOBAL, PASS, and ISOLNK. 


C.2. 3 IN?U 

TS 





• Calling 

Sequence: 

CALL RDDATA 
LAST) 

(ARRAY, 

TOP, IDATA, IDIM, 

Parameter 


Dimension 

In/Out 

Definition 

ARRAY 


TOP 


In 

Starting locution 
for field information. 

TOP 


1 


In 

Number of available 
words in ARRAY. 



Pareuneter 

Dimension 

In/Out 

Definition 

I DATA 

IDIM 

In 

Buffer for pixel 
value storage per 
scan line. 

IDIM 

1 

In 

Maximum dimension 




of IDATA. 

LAST 

1 

Out 

Flag indicating 


(when set to a 
value of 1) that 
all classes for 
one set of control 
cards have been 
processed. 


C.2.4 OUTPUTS 
No additional outputs. 

C.2.5 STORAGE REQUIREMENTS 
Code-2316g Data-1350g 

C.2.6 DESCRIPTION 

Designated "other" and/or "unidentifiable" pixels are 
identified from field information. The radiance values are 
reset to values of 0 and/or 255. This serves as the flag in 
subsequent processing. 

The field description from cards are held internally. 

The trigger for DO/DU processing is added to labeled common PASS 
and set in RDDATA. Two special print symbols {* and #) for 
DO and DU pixels are provided. 


Special cluster numbers are assigned to DO and DU pixels. These 
numbers have values 1 greater and 2 greater than the number of 
regular clusters. 

C . 2 . 7 FLOWCHART 
N/A 


C.2.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




SOFTWARE SUBPROGRAM NO. 3 (SETUP/) 


C. 3 

The new capabilities of subprogram SETUP? are: 

1. read an optional DOTFILE control card for unit and file 
number, and set option trigger for subsequent use by 
ISODAT 

2. read an optional SUNANG control card, set the sun angle 
correction trigger for use by ISODAT, and extract the 
sun angles if those are present on the card. 

C.3.1 LINKAGES 

The subprogram SETUP? is called b^ the routine ISOCLS. Routine 

SETUP? calls card reader/interpreter routines FIND, NXTCIIR, 

NUMBER, and FLTNCL. 

C.3.2 INTERFACES 

Subprogram SETUP? interfaces with other routines through com- 
mon blocks GLOBAL, PASS, and ISOLNK. 

C . 3 . 3 INPUTr 


• Calling 

sequence: CALL SETUP? 

(ARRAY , 

TOP, ITIME) 

Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

TOP-CORBAS 

Out 

Stores means or a 
statistics Module 
STAT deck. 

TOP 

1 

In 

Maximum dimension 
of ARRAY (normally 
10,600) 


ITIME 


1 


In 


Count of number of 
times SETUP? called 


• New Control Cards; 


Keyword 

DOTFIL 


SUNANG 


SUNANG 


Parameter 

INPUT/UNIT=n , FILE=m 
(Default; Self- 
initialized start- 
ing) 


Function 

Defines the FORTRAN 
unit number n and file 
number m of the dot 
data file DOTFIL con- 
taining the starting 
vectors. 


TAPE 

(Default: No sun 
angle correction 
applied) 

^1 ' ^2 ^ ^ j 
n. are integers, 

(Default: No sun 
angle correction 
applied) 


Sun angles will be eX' 
tracted from the MSS 
data tape . 


^ n^ ] are the sun angles 
to be used in computing 
the sun angle correc- 
tions for use in the 
clustering algorithm. 

A sun angle must be in- 
put for each set of 4 
channels appearing on 
the CHANNEL control card. 


C.3.4 OUTPUTS 
No additional outputs. 


C.3.5 STORAGE REQUIREMENTS 
Code-2410g Data-671g 

C.3.6. DESCRIPTION 

The DOTFILE and SUNANG control cards are handled by ex- 
tending data vector INVEC to include the character strings 


DOTFIL and SUNANG. The unit and file number of DOTFILE are 
stored in GLOBAL labeled common. Sun angles appearing on the 
SUNANG card are placed in ISOLNK as sun angle card input. The 
sun angle correction triggers appear in ISOLNK and are set by 
SETUP?. 


C . 3 . 7 FLOWCHART 


N/A 


C.3.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




C.4 SOFTWARE SUBPROGRAM NO. 4 (SUNFAC) 


Subprogram SUNFAC is identical to subprogram SUNFCS. Please 
refer to the description of this subprogram in Appendix B, nxamber 
B-20. 

C . 4 . 1 LINKAGES 

Subprogram SUNFAC is called by routine ISODAT. 

C.4. 2 DESCRIPTION 

The sun angle gain corrections computed by SUNFAC are used in 
computing pixel/cluster mean distances in the clustering routine 
PSPLIT and in determining the largest standard deviations in 
routine ISODAT. 



C . 5 SOFTWARE SUBPROG'RAM NO. 5 (PSPLIT) 


Subprogram PSPLIT performs much of the calculative activity 
associated with the ISOCLS processor. For each iteration of the 
clustering algorithn, it assigns pixels to clusters, and computes 
the new mean and standard deviation vectors. 

C.5.1 LINKAGES 

PSPLIT is called by routine ISODAT. It calls the UNIVAC drum 
read/write routines in RANDIO. 

C.5.2 INTERFACES 

T’SPLIT interfaces with other routines through common blocks PASS 
and ISOLNK. 

C.5.3 INPUTS 


Calling sequence; 

PSPLIT (MEANS, 
AVP, AMN, MEN) 

STDEV, N, 

CLD, C, IPLACE, 

Parameter 

Dimension 

In/Out 

Definition 

MEANS 

( NOFEAT, MAXCLS) 

In/Out 

Cluster mean 
vectors 

STDEV 

( NOFEAT, MAXCLS) 

Out 

Cluster standard 
deviations 

N 

(MAXCLS) 

Out 

No. of pixels in 
each cluster 

CLD 

(MAXCLS, MAXCLS) 

In 

Inter-cluster 
distances, not 
used 

C 

(NOFEAT, NOPTS) 

Out 

Buffer holding 
pixel radiance 
values after drum 
read 




1 


I PLACE 

(NOPTS) 

Out 

Cluster no. pixel I 


NOPTS«buffer 


is assigned i 


size for number 
of pixels 


i 

AVP 

( NOFEAT, MAXCLS) 

Out 

Intermediate i 


NOFEAT=number 


results in calcu- j 


of channels 


lation of cluster ’ 

AMN 

(max = 30) 
MAXCLS=62 (max 
number of 
clusters) 


standard deviations • 

AMN 

(NOFEAT, MAXCLS) 

Out 

Intermediate 




results in calcu- 




lation of cluster 




means 

MEN 

(NOFEAT, MAXCLS) 

Out 

Not used ; in 
ISODAT call, the 
same storage as 
MEANS 


C . 5 . 4 OUTPUTS 
N/A 

C . 5 . 5 STORAGE REQUIREMENTS 
Code-1002g Data-137g 

C.5.6 DESCRIPTION 

As the pixel radiance vec _ors come in from drum, these are first 
tested for 0 or 255 values over all channels (DO/DU pixels) . DO 
and DU pixels are assigned cluster numbers greater than the 
ordinary cluster numbers. Any pixel not DO or DU is assigned to 
a cluster based on minimum distance from the means provided to 
PSPLIT. Sun angle correction can be applied to this distance 



computation by user option. After all pixels are assigned to 
clusters, the cluster means and standard deviations are recomputed. 

In order to save calculation time, PSPLIT has two basically identi- 
cal cluster assignment loops, one having incorporated the sun 
angle correction and the other without it. Thus, for runs without 
sun angle correction, this multiplication (with default values of 
1) is not performed. 

C.5.7 FLOWCHART 
N/A 

C.5.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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D.l SOFTWARE SUBPROGRAM NO. 1 (SETUP 4) 


The SETUP4 subprogram of the SELECT processor is modified to 
recognize the new option, "CLSWT", on the OPTION control card. 

If the new option is input, SETUP4 initializes the table of 
weights for interclass subclass pairs (=1.0) and for intraclass 
subclass pairs (=0.0). Also, SETUP4 is modified to include a new 
internal subprogram, subroutine INTWGT which accomplishes the 
function of initializing interclass subclass-pair weights. 

D.1.1 LINKAGES 

The SETUP4 subprogram is called from the processor driver pro- 
gram, SELECT. Subprogram SETUP4 continues to reference the 
following subprograms: EMFIL, BSTCHK, CRDSTA, EVLCHK, FIND, 

GRPSCN, NUMBER, NXTCHR, ORDER, PRTFLD, REDSAV, WGTCHK, and WGTSCN. 
An additional call to internal subroutine INTWGT is included in 
this modification of SETUP4 . (See section D.3 for description 
of subroutine INTWGT) . 

D . 1 . 2 INTERFACES 

SETUP4 interfaces with the driver program, SELECT, and other 
subprograms via the named common blocks FSL, GLOBAL, and INFORM. 

D.l. 3 INPUTS 

Calling sequence: CALL SETUP4 (ARRAY, TOP, STOPFG, JTIME, SUBRAY, 


Parameter 

SUESIZ) 

Dimension 

In/Out 

Definition 

ARRAY 

10600 

In 

Working array in 

TOP 

1 

In 

blank COMMON. 

The maximiam usable 


• 

Df-t" 

location in ARRAY. 
TOP=10600, set 
in MONTOR. 


STOPFG 


1 


Out 


Flag set*l to 
indicate that the 
$END* card has 
been read. 

JTIME 1 In Counter for the 

no. of times 
SETUP4 is called 
from SELECT. When 
JTIME=1, the line 
printer output is 
labeled with the 
HEADER array 
contents . For 
JTIMEt^I, the 
heading is not 
placed on the 
output. 

SUBRAY 12000 In Working array, 

provided from the 
SELECT program. 

SUBSIZ 1 In Maximum usable 

location in SUBRAY, 
set in SELECT 
program (=12000) . 


• Revised Control Card 

Keyword Parameter Function 

(col. 1) (col. 11-72) 

OPTION CLSWT SETUP4 sets a flag, WTKEY=1, 

to indicate that internal sub- 

program INTWGT is to initialize 
the table of weights in SUBRAY 
for interclass subclass pairs. 



D.1.4 OUTPUTS 


One addition to line printer output is provided by the modified 
SETUP4 subprogram. When the OPTION CLSWT control card is input, 
the printer listing of user requests will include: 

"USE AUTOMATIC INTERCLASS SUBCLASS WEIGHTS" 

There are no other changes to the output. The subclass-pair 
weights determined in SETUP4 are listed in the INTERCLASS 
SEPARABILITY TABLE , which is part of the usual printer output 
from the SELECT processor. 

D.1.5 STORAGE REQUIREMENTS 

Storage required: Data-791 Code-990 

D.1.6 DESCRIPTION 

When the OPTION control card is recognized by SETUP4 , a transfer 
is made to statement number 60, where there are tests for three 
possible options - STATS , PUNCH , or CLSWT . Unique internal 
flags are set for each of the options detached on an OPTION 
card. If CLSWT (actually, "C" as a first character) is deleted, 
the flag WTKEY is set = 1- The flag, WTKEY, is tested in order 
to output the appropriate listing of user requests (see section 
D.1.4, OUTPUTS), and to initiate a call to an internal (to SETUP4) 
subprogram, subroutine INTWGT. The call to INTWGT accomplishes 
the automatic initialization of interclass subclass-pair weights 
in a SETUP4 working array, SUBRAY. 

The flag, SETWGT, is the flag normally set by SETUP4 when a 
WEIGHTS card is read (SETWGT=2) . Both SETWGT and WTKEY are 
tested for the necessity of calling the subprogram, WGTCHK. If 
either SETWGT=2 or WTKEY=1 or both, subroutine WGTCHK is called. 




In the call to WGTCHK, the subclass-pair weights initialized by 
si:ibroutine INTWGT are passed in array SUBRAY and two unused storage 
locations, ARRAY (WGHS14) and ARRAY (WGHS14+1) , are used to pass 
the flags SETWGT and WTKEY to subroutine WGTCHK. Subroutine 
WGTCHK will place any input (via a WEIGHTS card) subclass-pair 
weights into the array, overriding weights placed there by sub- 
routine INTWGT, then perform its usual function of reducing the 
two dimensional subclass-pair weights array to a linear array and 
passing it to SETUP4 beginning at ARRAY (WGHS14) . At this point, 
control is returned to the processor driver program, SELECT, as 
usual . 


D.1.7 FLOWCHART 


N/A 


D.1.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 





D.2 software subprogram no. 2 .(WGTCHK) 


Subprogram WGTCHK of the SELECT processor, as modified, performs 
one of two functions: 

a. The usual capability of initializing all intersubclass pair 
weights =1.0, followed by overriding these preset weights 
with the weights in the input argximent, WGHT, for subclass 
pairs (and honoring the OTHERS weight assignment, for all 
unspecified subclass-pairs). The weight table returned to 
SETUP4 via the argument, WEIGHT, contains for each inter- 
subclass-pair either a weight =1.0, or a weight = the weight 
specified on a WEIGHTS card input to SETUP4, or a weight = 
the weight specified on an OTHERS card input to SETUP4. 

b. Bypass the initialization of subclass pair weights to 1.0 
and the OTHERS weight aooignment, if any. Store into the 
weight table returned to SETUP4 via the argument, WEIGHT, 
either the weight (s) for specific intersubclass-pairs, as 
named on WEIGHTS card(s) input to SETUP4, or the weights 
for interclass subclass pairs as initialized (=1.0) in 
internal subroutine INTWGT and sent to WGTCHK via the 
argviment, WRKRY. All remaining intersubclass-pair weights 
in the WEIGHT table will be set = 0.0 (via initialization 
in subroutine INTWGT) . 

The choice of functions to be performed by WGTCHK is derived from 
the value of a flag, WTKEY, sent to WGTCHK via calling argument 
WEIGHT (2) . 

D.2.1 LINKAGES 

Subprogram WGTCHK is called by the SETUP 4 subprogram. WGTCHK 
does not reference any lower level programs . 

D.2. 2 INTERFACES 

Subprogram WGTCHK does not use a common block. Interface with 
the SETUP4 subprogram is via the calling arguments . 




1 


D.2.3 INPUTS 
Calling sequence: 

Parameters 

WEIGHI 


CLSNAM 


NAMPR 


CALL WGTCHK (WEIGHT, Co^SNAM, NAMPR, WGHT, WPTR, 
WRKRY, NOCLS2) 


Dimension In/Out 

NOCLS2* In/Out 

(NOCLS2-D/2 

(NOCLS2=No. 

of subclasses) 


Def ■' nition 

Output table of 
weights for sub- 
class pairs. WEIGHT 
(1) is equivalent 
to ARRAY(WGHS14) 
in the calling 
program SETUP4 . 

On input, WEIGHT (1) 
is used to convey 
the flag, SETWGT, 
from SETUP4 to WGTCHK. 
WEIGHT (2) is used 


to convey the flag, 
WTKEY, from SETUP4 
to WGTCHK. 


N0CLS2 In 

(=no. of 
subclasses) 


Final set of sub- 
class names from 
SAVTAP file, after 
any grouping of 
subclasses performed 
due to GROUP control 
card. 


(2, WPTR) In Subclass names or 

name pairs input 
TO SETUP4 on 
WEIGHTS card(s). 
NAMPR (1, WPTR )= 
subnamo^ , 


EK-6 



WGHT 


WPTR 


WRKRY 


400 


In 


1 


In 


(NOCLS2, In 

IOCLS2) 

NOCLS2=No. 
of subclasses 


NAMPR(2,WPTR)= 

subname 2 i 
if input, other- 
wise=blank. 
NAMPR(1,WPTR) may 
also contain 
"OTHERS", from 
WEIGHTS OTHERS 
card. 

Subclass, sub- 
class pair, or 
"OTHERS" weights 
from WEIGHTS 
card(s) input to 
SETUP4 via the 
WGTSCN subprogram. 

The number of 
weights input to 
SETUP4 on "WEIGHTS" 
control card(s). 

Working array. If 
WTKEY=WEIGHT ( 2 ) +1 , 
contains initializa- 
tion weights for 
intersubclass pairs 
(=1.0), and interim 
storage of subclass 
pair weights, from 
WGHT array. If WTKEY 
(WEIGHT ( 2 ) ) =1 , WRKRY 
contains interclass 
subclass pair weights 
sent in from sub- 
routine INTWGT. 


N0CLS2 1 7n The total no. of 

subclasses for 
all classes. 

D . 2 . 4 OUTPUTS 

The modificaticns to WGTCHK makes no change in line printer 
or file output. 

D . 2 . 5 STORAGE REQUIREMENTS 

WGTCHK requires 296 locations in the Instruction Bank, 69 locations 
in the Data Bank, for a total of 365 locations. 

D.2.6 DESCRIPTION 

The calling sequence for WGTCHK is not modified. However, the 
use of the first argument in the calling sequence, WEIGHT, is 
modified so that on input to WGTCHK, WEIGHT (1) = SETWGT (from 
SETUP4) and WEIGHT (2) = WTKEY (from SETUP4) . Two flags internal 
to WGTCHK, by the same names (i.e., SETWGT and WTKEY), are 
initialised from WEIGHT (1) and WEIGHT (2) immediately upon entry 
to WGTCHK, 

If WTKEY+1, subprogram WGTCHK performs its usual function; i.e., 
pre-initializing in the working array, WRKRY, all intersubclass- 
pair weights either -=1.0, or set to the value on the WEIGHTS 
OTHERS card, if it is input to SETUP4 . Following pre-initial- 
ization in the working array, WRKRY, if SETWGT=2 the subprogram 
overrides the preset weights in the working array with weights 
from the input WGHT array, using subclass names in the input 
NAMPR array to direct the storage of input weights. When the 
working array is completed, the subclass-pair weight values from 
the working array are transferred to the weight table (WEIGHT) 
and output by WGTCHK via calling argument. 



If WTKEY=1, WGTCHK assumes that the working array, WRKRY, contains 
subclass-pair weights as initialized in subroutine INTWGT(SETUP4) . 
Pre-initialization of WRKRY is bypassed, and the flag, SETWGT, 
is tested. If SETWGT=2 , the weights in the working array are 
overridden by weights from the input WGHT array, using subclass 
names in the input NAMPR array to direct the storage of weights 
into the working array, WRKRY. If SETWGT?^2, the weights in WRKRY 
remain as initialized in subroutine INTWGT. The weights from 
WRKRY are then transferred to the weight table (WEIGHT) and out- 
put by WGTCHK via calling argument. 

D . 2 . 7 FLOWCHART 
N/A 


D.2.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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D . 3 SOFTWARE SUBPROGRAM NO. 3 (IHTWGT) 


The subprogram INTWGT is an internal subroutine in the SETUP4 
subprogram. The function of INTWGT is to set«1.0 the weights 
of all interclass subclass-pairs, and to set=0.0 the weights 
of all intraclass subclass-pairs. 

D.3.1 LINKAGES 

Subroutine INTWGT is called by the SETUP4 subprogram. INTWGT 
does not reference any lower-level subprograms. 

D.3.2 INTERFACES 

INTWGT interfaces with subroutine SETUP4 via the calling argu- 
ments. No common blocks are utilized by INTWGT. 

D . 3 . 3 INPUTS 


CALL INTWGT 

(SUBWGT, CLSUB, 

N0SUB2 , 

N0CLS2) 


Parameter 

Dimension 

Type 

In/Out 

Definition 

SUBWGT 

(N0SUB2, 

Real 

In 

Working array for 


NOSUB2) 



storage of inter- 
class subclass-pair 
weights and intra- 
class subclass-pair 
weights . Equivalent 
to SUBRAY (Wl) in 
SETUP4. 

CLSUB 

N0CLS2 

Integer 

Tn 

Input array con- 
taining the no. of 
subclasses in each 
class. Equivalent 
to ARRAY (N0SUB2) in 
SETUP4 . 


NOSUB 2 

1 

Integer 

In 

The total no. 
subclasses in 
classes 

of 

all 

NOCLS2 

1 

Intege : 

In 

The total no. 
classes 

of 


D.3.4 OUTPUTS 

All output of subroutine INTWGT is via the first calling argu- 
ment (See Section D.3.3 above.) 

D . 3 . 5 STORAGE REQUIREMENTS 

INTWGT is internal to subroutine SETUP4. Approximately 50 storage 
locations are required by subroutine INTWGT alone. 

D.3.6 DESCRIPTION 

In SETUP4 when the flag, WTKEY=1, is set due to an input OPTION 
CLSWT control card, subroutine INTWGT is called by SETUP4 to 
perform the functions of setting interclass subclass-pair weights 
=1.0, and intraclass subclass pair weights = 0.0. 

INTWGT utilizes a portion of the working array, subray, to contain 
the subclass-pair weights returned to SETUP4. The portion of 
SUBRAY utilized by subroutine ];NTWGT is the same part of the array 
that is utilised as a working array (WRKRY) in subroutine WGTCKK. 
Within INTWGT, the array is known as SU3WGT, dimensioned (NOSUB2, 
NOSUB2) , where NOSUB2 = the total number of subclasses, after any 
grouping required by input GROUP control card(s). 

The method of setting the subclass-pair weights in INTWGT is as 
follows : 

1. initialize the weight for all subclass pairs = 0.0 
SUBWGT (I,J) = 0.0 1=1, NOSUB2; J=1 , N0SUB2 

SUBWGT (J,I) = 0.0 


2 


total number of classes 


. given that N0CLS2 = 
and ^ CLSUB(Kl) 

and = CLSUB(K2) 

then for K1 = 1, 2 , 
and K2 = Kl+1, 

the subclass numbers 
Kl-1 


= number of subclasses in class 
= number of subclasses in class 

3, . . N0CLS2-1 
Kl+2, . . N0CLS2 

, I, for class Kl are as follows 


K1 

K2 


i 


T — n 



1, 2, 3, 


N 


Kl 


(N = 0) 


likewise, the subclass numbers 
K2-1 

Jj = j + 3 = If 2, 3, . 

M=1 


J, for class K2 are as follows: 


f 


N 


K2 


= Number subclasses in class 1) 


3. SUBWGT (I,J) = 1.0 
SUBWGT (J,I) =1.0 


4. Return to SETUP4, with interclass subclass-pairs = 1.0 and 
intraclass subclass-pairs = 0.0, in SUBWGT. 


D . 3 . 7 FLOWCHART 
N/A 

D.3.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 
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E.l SOFTWARE SUBPROGRAM NO. 1 (REDIF2 ) 

Subprogram REDIF2 of the CLASSIFY procossor wai-- 
to, (1) accept a revised APRIORI control card and (2) set 
a flag (APRFLG) to indicate that in the subprogram SETUP2, the 
apriori probability values for each subclass are to be com- 
puted, using subclass population data from the statistics 
file, SAVTAP. 

E. 1.1 LINKAGES 

Subprogram REDIF2 is called by subprogram SETUP2. Subprogram 
REDIF2 calls the following subprograms: 

BMFIL, CATSCN, CRDSTA, FIND, FLTNUM, GRPSCN, NUMBER, NXTCHR, 
ORDER, and TDATE. 

F . 1 . 2 INTERFACES 

The common blocks used by subprogram REDIF2 are: CLASS, GLOBAL, 

INFORM. The variable, APRFLG, in common block CLASS was 
involved in the modification to subprogram REDIF2 . 

E.l. 3 INPUTS 

• Callixig sequence: Call REDIF2 (ARRAY , TOP, APRIOR, KATNO, 

BMATRX, PRIORI) 

Parameters Dimension In/out Definition 

ARRAY 10 600 In Variably dimensioned 

working storage 

TOP 1 In Maximum usable location 

in ARRAY - TOP = 10, 600 

In Storage for a priori 

probability values — 
Initialize = 0.0 in 
REDIF2 


APRIOR 


60 


Parameters Dimension In/out 


Definition 


KATNO 

BMATRX 


PRIORI 


60 In 

(BMCOMB, Out 

BMFEAT) 

BMCOMB « no. 
of linear com- 
binations in 
B-raatrix 
BMFEAT = no. 
of channels 
used in 
B-matrix 

60 Out 


• Revised control card 
Keyword Parameter 


Unused in REDIF2. 

Storage location for 
B-matrix 


Storage for apriori 
probability values read 
from input APRIORI con- 
trol card. 


Function 


APRIORI FILE If the first nonblank character 

on the APRIORI card (parameter 
field) is an F- , the flag, 
APRFLG, is set = -777777. This 
value of APRFLG will initiate 
in SETUP2 the computation of 
subclass apriori probability 
values from the statistics file's 
subclass pixel population data. 


E . 1 . 4 OUTPUTS 

File output is unchanged by these modifications to REDIF2. 



Line printer output is modified to include the message “*** CLSFY/ 
REDIF2 - BAD INPUT ON APRIORI CARD - DEFAULT APRIORI PROBABILITY 
VALUES WILL BE USED. This message is printed immediately follow- 
ing the listing of the APRIORI control card as is usually done, 
and is only printed when the first character on the APRIORI con- 
trol card is neither a legitimate "F", nor a number in the range 
of 0, 1, . . . , 9. 

E. 1 . 5 STORAGE REQUIREMENTS 

Subroutine REDIF2 requires 708 locations in the Instruction Bank, 
343 locations in the Data Bank, for a total of 1051 locations. 

E.1.6 DESCRIPTION 

REDIF2 modifications include; (1) a test for the character "F" 
as the first non-blank character on t.ie APRIORI control card, 

(2) the initialization of the variable, APRFLG = -777777 if the 
character "F" is encountered. 

If the first nonblank character is not an F , REDIF2 proceeds 
to the usual reading of apriori probability values from the 
parameter field of the APRIORI card. If the first character is 
not a legitimate numeric (number) , an error message is printed 
(see section E.1.4). In the case of an error on the APRIORI 
control card the flag variable, APRKEY is set = 0. 

E . 1 . 7 FLOWCHART 
N/A 


E.1.8 LISTING 

The listing is available in the Data Techniques Laboratory, JSC 
Building 17. 


E.2 SOFTWARE SUBPROGRAM NO. 2 (SETUP2) 


Subprogram SETUP2 tests the value of the flag, APRFLG, which 
is returned from subprogram REDIF2 via common block CLASS and 
if APRFLG ® -777777, SETUP2 uses subclass pixel population 
data obtained from the statistics file, SAVTAP, to compute 
the subclass apriori probability values which are stored in the 
APRIOR array, and subsequently used in classification. 

E.2.1 LINKAGES 

Subprogram SETUP2 is called by the CLASSIFY processor driver 
program, CLSFY. Subprogram SETUP2 calls the following sub- 
programs: FIND, NUMBER, REDIF2, REDSAV, and WRTBM. 

E.2. 2 INTERFACES 

Subprogram SETUP2 uses the named common blocks CLASS, 

GLOBAL, and INFORM. The variable, APRFLG, in common block 
CLASS was involved in the modifications to SETUP2. 

E.2. 3 INPUTS 

• Calling sequence; CALL SETUP2 ( ARRAY, TOP, FLDFLG, APRIORI , 

BMATRX , KATNO) 


Parameter 

Dimension 

In/out 

Definition 

ARRAY 

10/ 600 

In 

Variably dimensioned 
working storage. 

TOP 

1 

In 

Maximum usaole location 
in ARRAY - set = 10/ 600 

FLDFLG 

1 

In 

Unused in SETUP2 

APRIOR 

60 

Out 

Storage for the subclass 


apriori probability 
values . 



Parameter 


In/out 


Definition 


BMATRX 


KATiTO 


Dimension 


(BMCOMB, OUT 

BMFEAT) 

BMCOMB = no. 
of linear com- 
binations in 
B -matrix 
BMFEAT = no. 
of channels 
used in the 
B -matrix. 


60 


Out 


Storage for the B - 
transformation matrix 


Class-category 
correspondence . 


E . 2 . 4 OUTPUTS 

There were no changes in either file output or line printer 
output as a result of these modifications to SETUP2. 


E.2.5 STORAGE REQUIREMENTS 

Storage required: Data-looOo Code- 1442 ^ 

8 8 


E.2.6 DESCRIPTION 

The modifications to SETUP2 include: 

1. A test on the flag, APRFLG = -777777, to trigger a 
transfer to the coding for computation of subclass 
apriori probability values using subclass pixel popu- 
lations available in the array, KEPPTS (common block 
INFORM) . 

2. The computation and storage of subclass apriori proba- 
bility values, in the following manner: 

Total no. of subclass pixels 

NOSUB2 

TKEPTS = ^ KEPPTS (I) 

1=1 







apriori probability value for subclass I 


APRIOR(I) 


KEPPTS(I) 

TKEPTS 


I = 1,N0SUB2 (N0SUB2 = total no. of subclasses) 


E.2.7 FLOW CHART 
N/A 


E.2.8 LISTING 

The listing is available in the Data Techniques Laboratory, 
JSC Building 17. 



APPENDIX F 


SUBPROGRAM DESCRIPTIONS 




APPENDIX F 


TABLE OF CONTENTS 


F.l Software 
F.2 Software 
F.3 Software 
F.4 Software 
F.5 Software 
F.6 Software 
F.7 Software 
F . 8 Software 


for Subprogram 
for Subprogram 
for Subprogram 
for Subprogram 
for Subprogram 
for Subprogram 
for Subprogram 
for Subprogram 


No. 1 (SETUPS) 
No. 2 (DTCHK) 
No. 3 (REDIF3) 
No. 4 (DSPLY2) 
No. 5 (DOTPCT) 
No. 6 (DOTSUM) 
No. 7 (DSPLAY) 
No. 8 (TAP LAB) 


F.l SOFTWARE SUbPROGKAM NO. 1 (SETUP3) 


The modifications to SETUPS adapt SETUPS to LACIE Procedure 1 
image processing requirements. The modifications accomplish 
the following purposes: 

1. Implement the input of dot data parameters from the DOTFIL 
file created by the DOTDATA processor, when the new DOTFIL 
control card is input. 

2. Utilize existing storage in ARRAY normally occupied by train- 
ing field information to store the dot data parameters 

3. Perform a check for agreement between the set of dot category 
names and the set of classification category names from the 
input classification file, MAPTAP, and 

4. Order the dot data in storage by ascending line number with 
the group of dots on a given line ordered by ascending sample 
number . 

An added dot performance flag, DOTKEY and an associated error 
flag, DOTERR, IS TESTED. If DOTKEY>0 and DOTERR=0, activi- 
ties related to dot performance summary processing is 
initiated in SETUPS. If either DOTKEY=0 or DOTERR>0, dot per- 
formance summary processing is not initiated. 

If dot data processing is initiated and if there is any mismatch 
of category names between analyst-labeled dot categories (DOTFIL) 
and classification categories (MAPTAP) , dot data processing is 
be terminated, a message written to warn the user (see F.l. 4 (2) 
below) and the usual DISPLAY processor functions performed as 
detailed by the input control cards. Also, other error checks 
are made on the input DOTFIL data, before proceeding with dot 
processing (see items 5 and 7 in section F.l. 4). 


F.1.1 LINKAGES 


Subprogreun SETUPS is called by the DISPLAY processor driver pro- 
gram, DSPLAY. Subprogram SETUPS makes calls to the following 
subprograms: FIND, NUMBER, REDIFS, WRTFLD, FSBSFL, internal 

subprograms SETKEY and FIXFLD and three new subprograms TAPLAB, 

RDDOTS, and internal sxabprogram DTCHK. 

F . 1 . 2 INTERFACES 

Subprogram SETUPS references two common blocks, DISPL and GLOBAL. 

The modifications to SETUPS resulted in the addition of two 
variables to the common block, DISPL. The two added variables 
are DOTKEY and DOTERR. 

F.l.S INPUTS 

Calling sequence: CALL SETUPS (ARRAY, TOP) 

Parameter Dimension In/Out Definition 

ARRAY 10 600 In Variably dimen- ; 

sioned working S 

array . I 

j 

TOP 1 In Maximum usable 

location in ARRAY. 

Set-10 600 in the 
EOD-LARSYS execu- | 

i 

txve program, | 

MONTOR. I 

1 

I 

F.1.4 OUTPUTS i 

Modifications to SETUPS include eight additions to line printer i 

output : 

1. In the user input summary of requests the following message J 

was added: 



PRINT DOT DATA PERFORMANCE SUMMARIES FOR DOT DATA FROM 
FORTRAN UNIT XX, NO. XX AND TAPE (OR FILE) CCCCCCtll 

If OPTION OUTLINE control card is input, the following addi- 
tional line is printed: 

OUTLINE THE DOTS ON THE CLASSIFICATION MAP 

2. If each category riame from the input dor data file, DOTFIL, 
does not agree with one of the category names from the input 
classification results file, MAPTAP, the following message 
will be printed (from internal subprogram DTCHK) : 

DOT DATA PERFORMANCE SUMMARIES WILL NOT BE PRODUCED - THE 
CATEGORY NAMES FROM MAPTAP AND DOTFILE DO NOT MATCH - 
CATEGORY NAMES FROM MAPTAP ARE : MCATl , MCAT2 , MCAT3 , ... 

CATEGORY NAMES FROM DOTFILE ARE: DCATl , DCAT2 , DCAT3 , ... 

3. If each cateogry name from DOTFIL matches each category name 
on MAPTAP, or if each category name on DOTFIL matches one of 
the categories on MAPTAP and there are additionax categories 
on the MAPTAP which do not match a DOTFIL category name, the 
following message is printed (from internal subprogram DTCHK) : 

THERE IS OR IS NOT 1-1 CORRESPONDENCE BETWEEN THE CATEGORY 
NAMES FROM DOTFILE - DCATl , DCAT2 , ... 

AND THE CATEGORY NAMES FROM MAPTAP - MCATl , MCAT2 , MCAT3 , ... 

4. If the number of categories (NOCAT) on the input classification 
results file, MAPTAP, is zero(=0) (i.e., category classifica- 
tion not performed) , the following message is printed: 

*** CLASSIFICATION BY CATEGORY (ON MAPTAP) IS REQUIRED IN 
ORDER TO PROCESS THE DOT DATA *** 

*** DOT PERFORMANCE SUMMARIES WILL NOT BE OUTPUT *** 



5. If the number of categories (NDCAT) on the input dot data 
file, DOTFIL, is zero (*0) , the following message is printed: 

*** SETUP3 ** FROM DOTFILE, THE NO. OF DOT CATEGORIES-XX 
**** dot PERFORMANCE SUMMARIES WILL NOT BE PROVIDED **** 

6. If the total working storage available in ^RRAY (10600 
locations) is exceeded by the usage of ARRAY in SETUPS, the 
amount of excess storage locations attemted is output in the 
message : 

***** DISPLAY/SETUP3 - CORE OVERFLOW (TOP-TOP2) BY 
XXXX — EXECUTION TERMINATED ****** 


7. If the number of dots on the input dot data file, DOTFIL, 
exceeds 250 dots, the following message is printed: 

**** DSPLAY/DSPLY2 **** NO. OF DOTS XXXX EXCEEDS THE 

MAX. ALLOWABLE (250) 

**** DOT PERFORMANCE SUMMARIES WILL NOT BE PRODUCED **** 

8. After reading the dot data from DOTFIL, vising subroutine 
RDDOTS to read the file, a printout of the dot file contents 
is provided, as follows: 

DOT FILE INFORMATION: 

NO. DOT CATEGORIES = XX 
TOTAL NO. DOTS = YYY 
DOT CATEGORY NAMES: 

DCATl , DCAT2 , DCAT3 , ... 


DOT RECORD: 
DOT NO. 

SAMPLE 

LINE 

TYPE 

CATEGORY 

1 

XXA 

YYA 

2 

1 

2 

XXB 

YYB 

2 

2 

3 

XXC 

YYC 

2 

2 

• 

• 

• 

• 

• 

. 

• 

• 

• 

• 


• 

• 

. 

. 


m 


F.1.5 STORAGE REQUIREMENTS 


SETUPS requires 1530 locations in the Instruction Bank, 749 
locations in the Data Bank, for a total of 2279 locations. 

F.1.6 DESCRIPTION 

The modifications to subprogram SETUPS are as follows: 

1. The storage allocation in ARRAY for training fields {TRNSAV, 
TRNFLD, TRNVER) is used for final storage of the dot data 
information input by the subprogram RDDOTS. This storage 
area is used after the usual printout of SAVED TRAINING 
FIELDS. Treating each dot as a one pixel training field, 

the base addresses FLDSV2 , VERTX2 , and FIELD2, for three areas 
in ARRAY are recomputed, based on the total number of dots 
(<250) read from DOTFIL. 

Base addresses for ARRAY storage are computed as follows: 

FLDSV2 = TOPI, where TOPI is the first available location 
in ARRAY, above the area in ARRAY used for test 
field information. 

FIELD2 = FLDSV2 + 1000 (allowance for information on 250 dots 
in TRNSAV (4,250) 

VERTX2 = FIELDS + 1250 (allowance for information on 250 
dots in TRNFLD (5,250)) 

TOPS = VERTX2 + 500 (allowance for information on 250 dots 
in TRNVER (2,250)) 

2. The tests on TRNKEY are modified to include a test on the 
flag, DOTKEY (dot data processing flag) . 

3. Initiate the output of a correct list of user-requested 
options. (See section F.1.4(l)). 



The dot data perfomance summary flag, DOTKEY=l, is used to 
initiate input of dot data from the DOTFIL, and also a com- 
parison of category names input from the classification 
results file, MAPTAP, against the category names read from 
the dot data file, DOTFI*.,. If there is a one-for-one match 
on category names, dot data performance summaries will be 
produced by DISPLAY. If there is a mismatch between the two 
sets of category names, the message described in section 
F.1.4 (2) above will be written, the DOTKEY flag reset = 0, 
TRNKEY set = 0, and the error flag DOTERR set = 0. Dot data 
processing will not be done by the DISPLAY processor, when 
DOTKEY = 0 or DOTERR > 0. 

TOP2 is the beginning address in ARRAY for a scratch area, 
which is used to temporarily contain the input dot data from 
DOTFIL by subroutine RDDOTS . The maximum size required is 
1000 locations, for DOTS (4, NDOTS) . Subroutine RDDOTS 
initializes DOTS (1,1) = dot sample number, DOTS (2,1) = dot 
line number, DOTS (3,1) = dot type, and DOTS (4,1) = dot 
category number (on DOTFIL) . Subroutine RDDOTS also initia- 
lizes the total number of dots (TOTVT2) , the number of dot 
categories (NDCAT) , and the dot category names (DOTCAT) . 

The existing array, CARD (62) , in SETUPS is used to temporarily 
contain the dot category names read from DOTFIL by subroutine 
RDDOTS (i.e., DOTCAT, is equivalenced to CARD). 

An internal (to SETUPS) subprogram, subroutine DTCHK, checks 
the dot category names from DOTFIL against the classification 
category names from MAPTAP. If each dot category name is 
present in the set of classification category names, DTCHK 
orders the dot data in DOTS (4, NDOTS) (i.e., ARRAY (TOP2) ) 
such that the dots are in sequential ascending line number 
order, with all dots on a given line ordered by ascending 


pixel numbers. After ordering i-t-j dot data, DTCHK stores the 
dot data into the areas in ARRAY usually reserved for training 
field data. The training field data input from MAPTAP is 
overwritten in the process. 


The training field storage areas in ARRAY are initialized 
with dot data information as follows: 


a. Storage is allocated such that 

TRNSAVd,!) = ARRAY (FLDSV2) , 

TRNFLD(1,1) ® ARRAY (FIELD2) , 

TRNVER(1,1) = ARRAY (VERTX2) 

b. The training field-related parameters, T0TVT2 and N0FLD2, 
are set = the number of dots 


c . Dot data i j moved from the temporary storage in DOTS 
(4,ND0TS) to the training field storage areas in ARRAY 
in the following manner: 


TRNS.AV(1,I) 

TRNSAV(2,I) 

TRMSAV(3,I) 

TRNSAV (4,1) 
TRNVER(1,I) 
TRNVER (2,1) 

For 1=1,2 


= unused in SETUPS 

= MAPTAP category number corresponding to dot 
I analyst-labeled category name 
= dot type (usually = 2, for bias-correction 
dots) 

= 1 = number of "vertices" 

= dot sample nvimber 
= dot line number 

3,..., N0FLD2 (= total mamber of dots). 


Internal subprogram FIXFLD initializes TRNFLD as usual 
to contain the coordinates of the rectangular area 
surrounding each dot, such that 


TRNFLD (1,1) 
TRNFLD (2, I) 
TRNFLD (3, I) 


starting line number 
ending line number 
starting sample number 



TRNFLD(4,I) - ending sample number 

TRNFJ.i‘{5,I) ■ pointer to the location in TRNVER for the 
vertex of the field (dot) 

TRNFLD is later utilized to outline the dots on the 
classification map. 

8. The size of the performance table PCTAB is set = N0FLD2 * 
total number of dots (i.e., PCTSZ « N0FLD2 if DOTKEY > 0) and 
PCTAB utilizes the scratch area in ARRAY, beginning at ARRAY 
(TOP2) , as usual. 

9. After the ordering of the dot data, initialization of TRNSAV, 
TRNVER, and TRNFLD, the usual call to the subprogram SETKEY 
(if DOTKEY > 0) , will be bypassed, and SETUPS will set FLDKEY 
0, DOTKEY = mrniber of dot categories, and will take the 
normal return to DSPLAY. 

F . 1 . 7 FLOWCHART 


N/A 


F.1.8 LISTING 

The listing for SETUPS is available in the Data Techniques 
Laboratory, JSC Building 17. 




F.2 SOFTWARE SUBPROGRAM NO. 2 (DTCHK) 


DTCHK IS a new internal subroutine in subroutine SETUPS. The 
functions of this routine are; 

1. To check the correspondence between DOTFILE's and MAPTAP's 
category names. 

2. To order the dots by ascending line number, and for all dots 
on a given line, by ascending sample number. 

F.2.1 LINKAGES 

The subprogram, DTCHK, is called by SETUPS and is internal to 
SETUPS. The subprogram is not available to be referenced outside 
of SETUPS. 

F.2. 2 INTERFACES 

The interface with SETUPS is by means of three calling arguments 
and five global variables. 

F.2, 3 INPUTS 

Calling sequence: CALL DTCHK (DOT, TRNSAV, TRNVER, DCATNM, 

NDOTS, NDCAT) 




x'/ 




Parameter 

Dimension 

In/Out 

Definition 

DOT 

(4,ND0TS) 
where NDOTS * 
total number of 
dots (£ 250) 

In 

DOT (1,1) = sample 
no. 

(2.1) »line no. 

(3.1) *=dot type 

(4 . 1) =labeled dot 
(DOTFILE) cate- 
gory no. 

TRNSAV 

(4, NDOTS) 

Out 

TRNSAV (1,1) = 
unused 

(2.1) = category 
number 

(3.1) = dot type 

(4.1) = 1 = no. 


of "vertices" 


TRNVER 

(2,NOFLD2) 
where NOFLD2 = 
total no. of 
dots (£ 250) 

In 

TRNVER(1,I) = 
sample no. 
(2,I)=line no. 

DCATNM 

NDCAT 

(NDCAT=no. of 
dot categories) 

In 

The dot category 
names 

NDOTS 

1 

In 

Total no. of dots 
from DOTFILE 
(< 250) 

NDCAT 

1 

In 

Total no. of dot 
categories from 
DOTFILE (£60) 


F.2.4 OUTPUTS 


See ^he messages described in Section F.1.4 for the line printer 
outputs from DTCHK. 



F.2. 5 


STORAGE REQUIREMENTS 


DTCHK is internal to subroutine SETUPS. See Section F.1.5 for 
SETUPS storage requirements, including DTCHK. 

F.2.6 DESCRIPTION 

If there is a match in category names from MAPTAP for <.ac 
category name from DOTFIL , the type of correspondence (i.e., 

1-1 or NOT 1-1) is printed out along with each set of category 
names. The NOT 1-1 correspondence is for the case of more 
MAPTAP (class! f icat ion) categories than dot categories. 

If there is not a match for each dot category name in the set 
of classification (MAPTAP) categories, an error message is 
printed out. DOTKEY and TRNKEY is reset to zero (=0) and DOTERR 
is set to one (=1) . See Section F.1.4 for the message ..put. 

F.2.7 FLOWCHART 
N/A 


F.2.B LISTING 

The listing is available in the Data Techniques Laboratory, JSC 
Building 17. 


F-3 SOFTWARE SH lPRO fiR M,,.JlQ 3 (.RERIFll 

The modifications to subprogram REDIF3 result in the following 
additional functions being performed; 

1. A new control card, DOTFILE, is read and decoded 

2. A new flag, DOTKEY, is set and used in other subprograms to 
indicate that dot data classification performance summaries 
are requested. DOTKEY is added to the common block DISPL. 

3. A new flag, DOTERR, is initialized and used in other subpro- 
grams. DOTERR set to one if MAPTAP does not contain any dot 
category names. DOTERR is also added to COMMON/DI SPL/. 

F.3.1 LINKAGES 

Subprogram REDIF3 is called by subprogram SETUP3. Subprogram 
REDIF3 calls the following subprograms: CHIN, FIND, FLTNUM, 

LAREAD, NXTCHR. 

F.3.2 INTERFACES 

Subprogram REDIF3 references the common blocks DISPL and GLOBAL. 
Two new variables, DOTKEY and DOTERR are added to the common 
block DISPL. DOTKEY is tested in subprograms SETUPS, DSPLY2 , 
and DOTSUM, and used for initiating dot data processing. DOTERR 
may be reset c.r tested in subprograms SETUPS and DSPLAY. It is 
used to inhibit dot data processing, when an error condition has 
been detected. 

F.S.S INPUTS 

The calling sequence for subprogram REDIFS is unchanged by these 
modifications . 

Calling sequence: CALL REDIFS (TSTSAV, TSTFLD, TSTVER, 



VDTM) 


Parameter 


Dimension 


In/Out 


Definition 


TSTSAV (4,250) Out 


TSTFLD 


(5,200) 


Out 


TSTVER VDIM 


Out 


VDIM 


1 


In 


Used for test 
field or desig- 
nated field 
information: 

TSTSAV (1,N) = 
field name 
TSTSAV ( 2, N) = 
class number 
TSTSAV ( 3, N) = 
subclass 
TSTSAV ( 4, N) = 
no. of vertices 
for field N 
N = 1,2,3,... 
no. of fields 

Used for storage 
of coordinates 
of the rectangular 
area surrounding 
a test or desig- 
nated field. 

Used for storage 
of the sample, 
line coordinates 
of the vertices 
for all input 
test or designated 
fields . 

Total no. of 
locations allo- 
cated for storage 
of all test or 




or designated 
field vertices 
in the variably- 
dimensioned ARRAY 
in SETUPS. 

• New control card 
Keyword 
DOTFILE 

put of dot performance 
summaries in DSPLY2 . 

NOTE: If OUTPUT/UNI T=X, 
FILE=Y appears in the 
parameter field, it will 
be ignored. 


Parameter 


Function 


INPUT/UNlT=m, 

FILE=n 

(Default; UN1T=8, 
FILE=1) 


Implements the reading of 
dot data from the speci- 
fied (or default) unit and 
file in SETUPS and the out- 


F.S.4 OUTPUTS 

The line printer output lists the new control card at the time it 
is read by REDIFS. If MAPTAP does not contain any dot category 
names, on error message is printed, as follows; 

"*** CLASSIFICATION BY CATEGORY (IN MAPTAP) IS REQUIRED IN ORDER 
TO PROCESS THE DOT DATA ***" 

DOTERR is sot to one (=1) in this case. If INPUT is not found 
on the DOTFILE card, the default unit and file will be used to 
read the DOTFILE. If an error in format is sensed on the DOTFILE 
card following INPUT, the error message: 

"*** ERROR IN DOTFILE CARD-TERMINATE EXECUTION ***" 

is printed out and exit is called. 




F.3.5 STORAGE REQUIREMENTS 

REDIF3 requires 964 locations in the Instruction Bank, 411 

locations in the Data Bank, for a total of 1375 locations. 

F.3.6 DESCRIPTION 

The order of the modifications to subprogram REDTF3 are: 

1. The keyword array, OPT, is rcdimensioncd to 15 locations, 
and the new control card keyword, "DOTFIL", is added to OPT. 

2. The new dot data processing flag, DOTKEY = 0 is initialized, 
as well as TSTKEY (=^0), TRNKEY ( = 0), and DOTERR ( = 0). 

3. For match-up of keyword DOTFIL in OPT array versus input 
control card keyword, a transfer to a new statement number 
785 is added to the tiMnsfer list. 

4. Statement number 785 and additional coding to test for the 
number of MAPTAP categories (NOCAT) ^ 0 coding to set the 
new flag, DOTKEY = 1, and a transfer to statement number 789 
if NOCAT 0. At statement number 789 coding is added to 
set DOTERR = 1, TRNKEY 0. 

5. After existing statement number 870, if DOTKEY 0, or 
DOTERR '' 0, TSTKEY is set 0 (following the tost for input 
of test fields). Et^l lowing the test for DOTKEY '' 0, DOTERR ^ 
0, the usual test for "designated" field input is performed. 
They may be either "DESIGNATED OTHER" or "DESIGNATED 
UNIDENTIFIABLE." For dot data processing, the test field 
flag, TSTKEY, is turneii off (TSTKEY~0) . The designated field 
flan DESKEY is on (DESKEY=1), if and only if any fields read 
in are identified as "DESIGNATED" fields. 

F.3.7 FLOWCHART 

N/A 




F.3.8 LISTING 

The listing is available in the Data Techniques Laboratory, 
JSC Building 17. 


F<L6 


F-4 SOFTWARE SUBPROGRAM NO. 4 (nSPLV2) 

DSPLY2 is modified to test the dot performance flag, DOTKEY. 

If DOTKEY > 0, a call to a new internal subprogram DOTPCT is 
initiated to create the dot classification performance table, 
PCTAB. Also, a call to an additional new internal subprogram, 
DOTSUM produces the two dot data classification summaries. The 
capability to border each dot on the classification line printer 
map, as is normally done with training fields, is preserved. 

The dots are outlined on the classification map if OPTION 
OUTLINE control card has been input to the DISPLAY processor. 

F.4.1 LINKAGES 

Subprogram DSPLY is called by the DISPLAY processor driver pro- 
gram DSPLAY. Subprogram DSPLY2 calls the following subprograms: 
DESIG, FLDBOR, MAPHD, PCT, PRTSUM, WRTHED , WRTLN, PRTCOL, 

SETMRG (no longer effective), RWRITE, WRTHD, RREAD, and two 
new internal subprograms, DOTSUM and DOTPCT. 

F . 4 . 2 INTERFACES 

Subprogram DSPLY2 references two common blocks DISPL and GLOBAL. 
F.4.3 INPUTS 

The calling sequence for DSPLY2 is unchanged as a result of these 
modifications. The usage of the three training field information 
array (allocated in ARRAY) in the calling arguments - TRNSAV, 
TRNFLD, and TRNVER - are revised by these modifications. Also, 
the usage of the classification performance table PCTAB (also 
allocated to an area in ARRAY) , is modified for dot data perfor- 
mance tabulation. The revision to usage of the training field 
storage area is in the typo of information placed in the three 
areas, and is only active for the dot performance svimmary option 
(i.e., DOTKEY > 0). Otherwise, the training field storage areas 
are available for their normal usage. The main change to usage 




of PCTAB is the size (dimension) of the storage area used for 
dot classification tabulation. For data, PCTAB may be treated 
as a linear (one-dimensional) array, since there is only a one 
dot - one subclass correspondence possible in the classification 
of each dot. Once again, this usage of PCTAB is only activated 
by the request for dot performance summaries (DOTKEY > 0) , and 
normal usage of PCTAB for training or test fields remains 
available when not in the dot processing mode in the DISPLAY 
processor. 

Calling sequence: CALL DSPLY2 (TRNSAV, TRNFLD, TRNVER, TSTSAV, 

TSTFLD, TSTVER, PCTAB) 


Parameter 

Dimension 

In/Out 

Definition 

TRNSAV 

(4 ,NOFLD2) 

In 

Initialj.zed in 


NOFLD2 = total 


SETUPS as follows 


no. of dots 


if DOTKEY > 0) : 


TRNSAV (1,1) = 
unused (scratch) 
TRNSAV (2, I) = 

MAPTAP (classifi- 
cation) category 
no. , for the MAPTAP 
category matching 
the dot's category 
name . 

TRNSAV (3, I) = 

DOT type (=l,or 2) 
TRNSAV (4, I) = no. 
of "vertices" for 
DOT I, 

I = 1,2, ... , N0FLD2 


total no. of dots 


TRNPLD 


(5,NOFLD2) 


TRNVER 


TSTSAV 


In 


(2,TOTVT2) In 

(T0TVT2 = total 
no. of vertices 
for all dots = 
total number of 
dots 


(4,NOFLD3) In 

N0FLD2 = total 
no. of test or 
designated 
fields (^200) 


Coordinates of 
rectangular area 
surrounding each 
dot. Initialized 
in SETUP3 as follows 
TRNFLD(1,I) = 
starting line no. 
TRNFLD(2,I) = ending 
line no. 

TRNFLD(3,I) = 
starting sample no. 
TRNFLD(4,I) = enuing 
sample no . 
TRNFLD(5,I) = loca- 
tion in TRNVER for 
vertex of dot 1. 

I = 1,2, .. . ,N0FLD2 = 
total no. of dots 

Array storage for 
training field (dot) 
vertex coordinates: 
TRNVER (1,1) = sample 
niamber 

TRNVER (2, I) = line 
number 

I = 1,2,..., total 
no. of dots = NOFLD2 
TOTVT2 

Array storage for 
test or designated 
fields, equivalent 
to TRNSAV for 
training fields. 




TSTFLD 

(5,NOFLD3) 

In 

Used for test or 
designated fields, 
equivalent to TRNFLD 
for training fields. 

TSTVER 

(2,TOTVT3) 
(T0TVT3 = total 
no. of vertices 
for all test or 
designated 
fields 

In 

Used for test or 
designated fields, 
information equiva- 
lent to TRNVER for 
training fields. 

PCTAB 

(NGFLD2,1) 
N0FLD2 = total 
no. of dots 

In 

Storage area, 
beginning at ARRAY 
(PCTID3) , for dot 


performance tabula- 
tion (PCTAB(I,1) = 
classified subclass 
for dot I, 1 = 

1 , 2 , , N0FLD2 , 

upon return from 
internal subprogram, 
DOTPCT) . 

The input of DESIGNATED fields to DSPLY2, in the test field 
storage areas - TSTSAV, TSTFLD, and TSTVER - is preserved (not 
changed) , and if designated fields have been input to the 
DISPLAY processor (the flag DESKEY is set = 1) , the usual 
processing of the classification record to mark pixels in the 
designated area(s) is performed in DSPLY2, with the call to 
subroutine DESIGN, as usual. 

F . 4 . 4 OUTPUTS 

There are three additional messages printed from DSPLY2 , as a 
result of these modifications when in the dot processing mode 



EKDTKEY > 0. The dot classification performance table, PCTAB, is 
checked after initialization in internal subprogram DOTPCT. 

1. If the classification subclass number for any dot exceeds 

N0SUB3 (■ the threshold "subclass" number) , it is assumed 

that the dot has been classified "DESIGNATED", in which case 

the message is printed: 
sample, line 

"DOT (XXX, YYYFTS in the designated AREA" 

2. If the classification subclass number for a dot is zero (*0) , 
the dot has been found to have sample, line coordinates not 
included in the actual area classified, and the subclass for 
the dot set “ 0 in internal subprogram DOTACT. In this case, 
the message is printed: 

sample , line 

"DOT (XXX, YYY) IS NOT IN THE CLASSIFIED AREA" 

3. If the number of dot caLeg< ries matching MAPTAP (classifica- 
tion) categories is not equal to the total number of dot 
categories on the input DOTFILE, the message is printed (from 
internal subprogram DOTSUM) : 

"**** DSPLY2/DOTSUM — DISCREPANCY IN DOTFILE INFOR^IATION **** 
NO. OF DOT CATEGORY LABELS MATCHING MAPTAP CATEGORY NAMES = XX 
NO. OF DOT CATEGORIES IS GIVEN AS = YY" 

The classification map and the classification summary for the 
total area classified is output from DSPLY2, as usual. 

In addition, as a result of these modifications, following the 
usual classification map and associated classification summary, 
two dot classification performance siammaries are output, when in 
the dot processing mode (DOTKEY > 0) . The two dot performance 
summaries - one by dot category, the other by individual dot - 
are output from the new internal subprogram, DOTSUM. See Section 
F.6.4 for a detailed description of the output dot classification 
performance summaries. 


When in the dot processing mode {EK5TKEY > 0) , DSPLY2 ignores any 
"test" field processing - i.e., outlining of "test" fields on the 
classification map. 

When in the dot processing mode (DOTKEY > 0) and when the OPTION 
OUTLINE control card has been input to the DISPLAY processor, 
the flag, TRNKEY, will have been set » 1, in which case DSPLY2 
outlines each dot on the classification map with a border of 
askerisks {"*"), as is normally provided for training fields 
when not in the dot processing mode. 

F . 4 . 5 STORAGE REQUIREMENTS 

The subprogram, DSPLY2, including the two added internal subpro- 
grams, DOTPCT and DOTSUM, requires 2106 locations in the Instruc- 
tion Bank, 13340 locations in the Data Bank, for a total of 
15446 locations. 

F.4.6 DESCRIPTION 

The modifications to DSPLY2 are to add capabilities related to 
LACIE Procedure 1 image processing techniques. LACIE Procedure 
1 techniques involve dots, which are one-pixel fields, category- 
labeled by photo- interpreters. Normal (non-LACIE Procedure 1) 
capabilities and options are preserved in DSPLY2, when not in 
the dot processing (LACIE Procedure 1) mode. To achieve the 
additional capabilities, DSPLY2 is modified as follows: 

1. A test on the dot processing flags, DOTKEY and DOTERR is made, 
to bypass the usual initialization of the performance table, 
PCTAB, to zero. This is necessary due to the revised 
dimension of PCTAB when it is used for dot classification 
tabulations . 


F-2f2 




TRNSAV(4,I), I » N0FLD2 (no*, of dots) is set » 2 , to 
satisfy the requirement of the field bordering subprogreun, 
FLDBOR, for the number of vertices of a one-pixel field. 

The IR array, used in the classification performance to 
contain the classification record of three adjacent scan lines, 
is reused in the new internal subprogram, DOTSUM, as a scratch 
area for several arrays used internally in DOTSUM. The 
equivalence of these arrays to portions of the IR array is 
done in DSPLY2, thereby making these arrays global and unne- 
cessary to include in the calling sequence to DOTSUM. 

In the usual scan line-by-scan line classification performance 
loop in DSPLY2 , tests are added for dot processing (DOTKEY > 

0 and DOTERR ^ 1) , in order to go through additional coding 
which tests for the presence of dots on the current scan line, 
and provides a call to the new internal subprogram, DOTPCT, 
to build the dot classification tabulation in PCTAB. Tests 
are made before calling DOTPCT, to see that there are dots on 
the current scan line. 

The usual call to subprogram PCT, to build performance tables 
for either test or training fields, is by-passed when in the 
dot processing mode, but remains available as usual when not 
doing dot processing. 

The bordering of dots on the classification map is performed 
by subprogram FLDBOR, when the training field flag, TRNKEY, 
is set = 1. (This will be the case when the OPTION OUTLINE 
card is input to the DISPLAY processor) . Test field bordering 
is bypassed, by virtue of the test field flag, TSTKEY, having 
been set = 0, when in the dot processing mode. 

Following the call to subroutine PRTSUM, which prints the 
classification map, the dot classification table, PCTAB, is 
checked for dots which have subclass numbers greater than 
N0SUB3 (the thresholded subclass number) , or equal zero 
(=0) which is the case when internal subprogram DOTPCT 


determines that the dot coordinates are not compatible with 
the pixel coordinates of the classified area. Any dots in 
the classification table, PCTAB, having a subclass niamber 
greater than N0SUB3 is assumed to be either "DESIGNATED OTHER" 
(=NOSUB3+5) or "DESIGNATED UNIDENTIFIABLE" (=NOSUB3+4) , in 
which case the dot classification subclass is set = 0 in 
PCTAB. For each dot that is already set = 0 in DOTPCT, a 
message to the user is printed out to the effect that the 
dot is not contained in the set of pixel which were classified. 
For dots set = 0 due to being within a "designated" area, a 
message to the user is printed out identifying the dot as 
being in a designated area. These dots will not be considered 
in the dot classification summaries provided later in internal 
subprogram DOTSUM. 

8. Following the printout of the classification map by PRTSUM 

and the checking of dot classification described in (7) above, 
a call is made to the new internal subprogram, DOTSUM, which 
will provide the two dot classification summaries described 
in section F.4.4, using the dot classification tabulation in 
PCTAB as the basis for the summaries. 

F . 4 . 7 FLOWCHART 

The flowchart for the new ir ’:ernal subprogram, DOTSUM is provided 
in section F.6.7. The flowchart for the new subprogram, DOTPCT, 
is provided in section F.5.7. 

F.4.8 LISTING 

The listing is available in the Data Techniques Laboratory, JSC 
Building 17. 



F.5 SOFTWARE SUBPROGRAM NO. 5 (DOTPCT) 


DOTPCT is a new internal subprogram in subprogram DSPLY2. DOTPCT 
is called only if the DOTKEY flag is on (>0) . The functions of 
this subprogr 2 im are; 

1. If dot (sample, line) is not contained in the classified area, 
PCTAB (K,l) is set to zero. K refers to the appropriate 

dot (sample, line) . 

2. If dot (sample, line) is contained in the classified area, 
PCTAB (K,l) is set to its respective classified subclass 
number. 

F.5.1 LINKAGES 

Subprogram DOTPCT is called by DSPLY2 and does not reference any 
other subprograms, nor may be referenced by any other routine 
external to DSPLY2. 

P.5.2 INTERFACES 

The interface with subprogram DSPLY2 is by means of three calling 
arguments . 

F . 5 . 3 INPUTS 

Calling sequence: CALL DOTPCT (PCTAB, FIELD, IR) 


Parameter 

Dimension 

In/Out 

Definition 

PCTAB 

(N0FLD2,1) 

Out 

Table of subclass 


N0FLD2=no. of 


numbers for the dots. 


dots 


PCTAB(I,l)=classifi- 


cation subclass no. 
for dot I. 


FIELD 


(2,NOFLD2) 


In 


IR PTS 

PTS=no. of 
pixels on each 
classified scan 
line . 


Table of dot 
coordinates. FTKI.D 
(1,1) = pixel no. , 
FIELD(2,I)=line no. 
for dot I. (equiva- 
lent to TRNVER in 
DSPLY2) 

In The classification 

record for one scan 
line. IR(J)- subclass 
no. for pixel J. 


In addition to the calling arguments, subprogram DOTPCT utilizes 
global variables which have been initialized in DSPLY2 , or sent 
to DSPLY2 already initialized. This set of global variables 
required by DOTPCT are: 

1. PCTKEY = flag for initializing PCTAB to zero (=0) for all 

dots - this is done once on the first call to DOTPCT 
(i.e., first scan lino), and the flag PCTKEY is then 
reset in DSPLY2 (PCTKEY > 0) so that pre-initialization 
of PCTAB is bypassed after the first call to DOTPCT. 

2. SAMSTR = first pixel number in the classification record, IR. 

3. SAMINC = pixel increment, in the classification record IR. 

4. BCNT = starting position, in the dot "vertex" coordinate 

array, for dots having a lino number equal to the 
current classif icat ion scan lino number. 

5. ECNT - ending position in the dot "vertex" coordinate array, 

for dots havii.g a line number equal to the current 
classification scan line number. 


F.5.4 OUTPUTS 


The tabulation of subclass number for each dot is output by 
DOTPCT, via the calling argument, PCTAB. 

F . 5 . 5 STORAGE REQUIREMENTS 

Subroutine DOTPCT is internal to subroutine DSPLY2. The storage 
requirements for DOTPCT are included in the storage requirements 
for DSPLY2 (see section F.4.5). 

F.5.6 DESCRIPTION 

DOTPCT is called from DSPLY2 once for each classified scan line 
input from MAPTAP in DSPLY2 . Tests are made in DSPLY2 to 
determine that there are dots with the given scan line number, 
before calling DOTPCT. In DOTPCT, a check is made on the pixel 
number for each dot with the given scan line number, to ascertain 
that the dot is present in the actual pixels classified. This 
must be done due to the possibility of sample and/or line 
incrementation in the classified area by the CLASSIFY processor 
which generated the classification results file, MAPTAP. 

For those dots present in the given classification record (IR) , 
the subclass (number) at the dot's position in the classification 
record is extracted from the classification record (IR) and 
placed into the dot's position in the performance table, PCTAB. 

For those dots having pixel coordinates which are not in the 
classification record, the subclass number is placed in PCTAB 
as a zero (=0) . 

The loop for checking the dots and initializing the performance 
table is controlled by two global variables determined in DISPLY2 
BCNT and ECNT. BCNT is the beginning dot number for dots having 
a line coordinate equal to the current classified scan line 
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mamber, and ECNT is the ending dot number for dots having a line 
coordinate equal to the current classified scan line number. 

After FIELD (1,BCNT) through FIELD (1,ECNT) is processed (i.e., 
the dots BCNT through ECNT) a return is made to calling subprogram 
DSPLY2, with the classification performance for these dots 
tabulated and returned in the argument, PCTAB. 

F . 5 . 7 FLOWCHART 

See Figure F-5 for flowchart of subprogram DOTPCT. 


F.5.8 LISTING 

The listing of subprogram DOTPCT is contained within the listing 
of subprogram DSPLY2, which is available in the Data Techniques 
Laboratory, JSC Building 17. 
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F.6 SOFTWARE SUBPROGRAM NO. 6 (DQTSUM) 


Subprogram DOTSUM is a new internal subprogram added to the 
DSPLY2 subprogram in the DISPLAY processor. The functions 
performed by the subprogram are as follows: 

1. To compute analyst labeled category proportions over the 
classified area as follows: 

a. Uncorrected proportions for the analyst-labeled categories 

b. Bias corrected proportions for the analyst-labeled 
categories 

2. To compute the "alpha" table 

3. To output two classification performance tables as follows: 

a. Summary by dot (analyst-labeled) category showing, for 
each analyst-labeled category: 

1. ) Total number of dots labeled by analyst for the 

given category 

2. ) Percent of correct classification of analyst- 

labeled dots (i.e., ratio of number of dots labeled 
Cj, classified to the total number of dots 
labeled C^) 

3. ) Number of dots thresholded 

4 . ) Uncorrected category proportion 

5. ) Bias corrected category proportion 

6 . ) Alpha values 

7. ) Number of dots in the labeled category which were 

classified into each possible classification (MAPTAP) 
category. 

b. Summary by individual dot showing: 

1.) Coordinate of dot (sample, line) 



2. ) Analyst-labeled category name 

3. ) Classifier's identification of the dot (category 

and subclass name) 

F.6.1 LINKAGES 

Subprogram DOTSUM is internal to subprogram DSPLY2 . The subpro- 
gram does not reference any other subprograms, and may not be 
referenced by other routines external to DSPLY2 . 

F . 6 . 2 INTERFACES 

Certain variables from common block DISPL and DSPLY2 are utilized 
by subprogram DOTSUM, as global initializes. See the flowchart 
for DOTSUM (Figure F-6) for identification of these variables. 


F.6.3 INPUTS 
Calling sequence: 

Parameter 

SUBTOT 


THTOT 


PCTAB 


CALL DOTSUM (SUBTOT, THTOT, PCTAB, TRNSAV, 
TRNVER) 


Dimension 

NOSUB 3 
(no. of 

subclasses +1) 


NOSUB 3 


(N0FLD2,1) 
N0FLD2=total 
no. of dots 


In/Out 


In 


In 


In 


Definition 

Total no. of pixels 
classified into 
each of the sub- 
classes, and thres- 
holded total for all 
subclasses . 

Total no. of pixels 
thresholded in each 
subclass, over the 
total area classi- 
fied . 

For DOT performance, 
PCTAB (1,1) = sub- 
class for dot I 


TRNSAV 

(4,NOFLD2) 

In 

TRNSAV ( 1 , N) =unused 
TRNSAV ( 2, N) = labeled 
category number 
TRNSAV(3,N)==dot type 
TRNSAV ( 4, N) =2 (no. 
of vertices) 

TRNVER 

(2,TOTVT2) 

In 

Dot coordinates; 


T0TVT2=total 


TRNVER (1,1) = sample 


no. of dots 


coordinate, TRNVER 
(2,1) = line coordi- 
nate of dot I 


F.6.4 OUTPUTS 

The dot data classification svmunaries are provided in place of 
the usual test or training field summaries. Two dot classifica- 
tion summaries are output: 

• Dot category summary table 

• Individual dot summary table 

F . 6 . 5 STORAGE REQUIREMENTS 

The storage for DOTSUM is included in the total storage require- 
ments of subroutine DSPLY2 (see Section F.4.5). 

F.6.6 DESCRIPTION 

The subprogram functions as follows: 

1. Input via calling arguments and common block DISPL- 


a. SUBTOT (N0SUB3) = total number of pixels classified into 

subclass 1,2,..., NOSUB 2 , THRESH 

N0SUB2 = total number of subclasses 

N0SUB3 = total number of subclasses + 1 



b. THTOT (NOSUB3) = total pixels thresholded, by subclass 

c. PCTAB (NOFLD2,4) = the subclass no. for each dot i, i = 

1,2, . . . ,N0FLD2 = total number of dots 

d. TRNSAV (4,NOFLD2): 

TRNSAV(1,N) -- unused 

TRNSAV ( 2, N) * dot category number for the MAPTAP category 
matching the analyst-labeled category of 
the dot 

TRNSAV(3,N) == dot type ( = 1 or =2) 

TRNSAV ( 4, N) = 2 = number of "vertices" for dot 

e. From common block DISPL; NOCAT = nxamber of classification 
categories; CATNAM(61) = classification category names; 
SUBCAT (60) = subclass-to-category correspondence (SUBCAT 
(I) = M, means subclass number I belongs in category 
number M) ; NOFLD2 = total number of dots; NDCAT = DOTKEY = 
number of dot categories, from DOTFILE. 

Store dot type in TYPE (=TRNSAV(3,1) ) , initialize NDCAT = 
DOTKEY = number of dot categories and initialize NCLAT = 

NOCAT = number of classification categories, NCLCAT = NOCAT +1 . 

Move labeled category from TRNSAV (2, I) to TRNSAV (3, I), for 
1=1, NOFLD2 (=total number of dots) 


Transfer each dot "vertex" (sample, line coordinate) from 
TRNVER array to TRNSAV (1,1), TRNSAV (2,1) for 1=1, 
NOFLD2 dots [TRNSAV (1,1) = dot sample, TRNSAV (2, I) = dot 
line number, TRNSAV (3, 1) = labeled category nuitjber, 1 = 1, 
NOFLD2 = total number of dots] . 


Compute total number of pixels classified in the total area 
classified, except for thresholded pixels and pixels in the 
"designated" area: NOSUB2 


PIXTOT 


SUBTOT(i) 


Compute the total number of pixels classified into each 
analyst-labeled category; 

for I » 1,2,..., N0SUB2 (total number of subclasses), 
LBLTOT(K) - LBLTOT(K) + SUBTOT(I), 
where K “ SUBCAT (I) » category 
corresponding to subclass I. 

Compute dot totals classified, by classification category: 
Initialize CLTOT(K) * 0, K * 1 , 2 , . . . ,NCLAT (=NOCAT) . Then, 
for I » 1, NOFLD2, SUBCL ^ (PCTAB(I,1), 1 < SUBCL < NOSUB2 , 
CAT = SUBCAT (SUBCL) ; CLTOT(CAT) * CLTOT(CAT) + 1 

Based on occurrence in the dot data array, TRNSAV, initialize 
an ordered set of dot category labels DOTLBL(J,l) = 
CATNAM(LBL) = category name and D0TLBL(J,2) = LBL = category 
number of labeled dot for J = 1, 2 , . . . ,NDCAT, and LBL = 
TRNSAV(3,I), I = 1,2,..., no. dots, such that DOTLBL 
(J+1,2) > DOTLBL(J,2). 

Compute the uncorrected proportions of dot categories over 
the total area classifeid; PUNC(I) = LBLTOT ( J) /PIXTOT^ 100 ; 
1=1, NDCAT; J = DOTLBL(I,2). 

Initialize TRNSAV(4,I), 1=1, NOFLD2 (number of dots) to 
contain the classification category (number) for dot I: 

SUBCL = PCTAB(I,1), I = 1 , 2 , . . . ,NOFLD2 ; if SUBCL = NOSUB3, 

CAT = NOCLCAT (=NOCAT + 1); if SUBCL NOSUB3, CAT = SUBCAT 
(SUBCL); TRNSAV(4,I) = CAT (classified category number). 

Initialize LABCLS(I,J) = 0, I = 1 , 2 , . . . ,NODCAT; J = 1,2,..., 
NCLCAT (=NOCAT + 1) 

Form a table showing (a) the labeled dot category, and (b) 
the total number of dots with a given label that have been 
classified into each classification category: 

Initialize LABCLS(I,J) =0, 1=1, NDCAT, J = 1, NCLCAT; then 
for J = 1, NOFLD2, LB = TRNSAV (3,J), CLS = TRNSAV (4,J): 
LABCLS (LB, CLS) = LABCLS (LB, CLS) + 1. 


13 . 


Compute alpha table for labeled versus classified dots; 
alpha (I,J) » LABCLS (J,I)/CLTOT(I) *100, I » 1, NOCAT classi- 
fication categories, J * 1, NDCAT labeled categories; set 
alpha (I,J) « 0, where I > NDCAT. 

14. Compute corrected proportions for labeled dot categories: 
for 1*1, NDCAT (no. of dot categories) , 

NDCAT 

PCORR(I) * z: [alpha(J,D * PUNC(J)] * 100 

J*1 

15. Tabulate the total number of analyst-labeled dots in each dot 
category: 

initialize TOTLBL = 0; then for 1=1, NOFLD2 = total number 
of dots and LBL = TRNSAV(3,I) and for J where DOTLBL(J,2) = 

LBL; TOTLBL (J) = TOTLBL (J) + 1. 

16. Tabulate the number of dots in each category which were 
thresholded; initialize TOTTHR(J) =0, 1=1, NDCAT; then 
for 1=1, NDCAT * no. of dot categories, TOTTHR(I) = 

LABCLC (I,NCLCAT) 

17. Print first dot classification performance summary, by 

category: for 1=1, NDCAT, 

a. labeled dot category name = DOTLBL(I,l) 

b. total number of analyst-labeled dots in each category = 
TOTLBL ( I ) 

c. percent correct classification = LABCLC (I , I) /TOTLBL (I) *100 . 

d. number of dots thresholded in the analyst-labeled 
category = TOTTHR(I) 

e. uncorrected labeled category proportion = PUNC(I) 

f. bias-corrected labeled category proportion = PCORR(I) 

g. classification summary and Ujj for labeled category I: 
for J = 1, NCLCAT; LABCLS(I,J), ALPHA(J,I) 


18. Set SUBNAM(N0SUB3) - 'THRESH' 

19. Print second dot classification performance summary, by 

individual dot: for I ■ 1, N0FLD2 ■ total number of dots. 

a. dot identification “ (sample, line) * TRNSAV(1,I), 
TRNSAV(2,I) 

b. labeled category name » CATNAM[TRNSAV(3, I) ] 

c. classified category name = CATNAM[TRNSAV(4,I) ] 

d. classified subclass = SUBNAM [ PCTAB { I , 1 ) ] 

e. If dot is thresholded, in designated area, or not in the 
area classified, substitute "******•’ for both category 
and subclass neunes to be output for the dot. 

See figure F-6 for flowchart of internal subprogram DOTSUM. 


F . 6 . 7 FLOWCHART 
N/A 


F.6.8 LISTING 

The listing for DOTSUM is contained in the listing of subprogram 
DSPLY2, which is available in the Data Techniques Laboratory, 

JSC Building 17. 
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F.7 SOFTWARE SUBPROGRAM NO. 7 (DSPLAY) 


DSPLAY is the driver program for the DISPLAY processor. Minor 
modifications were required to incorporate LACIE Procedure 1 image 
data processing options. The primary impact of the modifications 
to DSPLAY are described in Section F.7. 4 and F.7. 6 below. 

F.7.1 LINKAGES 

DSPLAY is called by the EOD-LARSYS executive program, MONTOR. 
DSPLAY calls SETUPS, DSPLYl , EMTHRS , DSPLY2 , and PRTPCT. 

F.7. 2 INTERFACES 

The display program references the DISPL common block. 

F.7. 3 INPUTS 

Calling Sequence: CALL DSPLAY (ARRAY, TOP) 


Parameter 

Dimension 

In/Out 

Definition 

ARRAY 

10,600 

In 

Variably dimensioned 
working array 

TOP 

1 

In 

Max. usable location 
in array. SET~10,600 
in the EOD-LARSYS 
executive program, 
MONTOR . 


F.7. 4 OUTPUTS 

VJhen not in the dot (LACIE Procedure 1) processing mode, a call 
is made to EMTHRS from DSPLAY, when tost or training fields have 
been input to the DISPLAY processor. EMTHRS calculates the 
emperical thresholds, using data from defined test or training 
field storage areas in the largo working array, ARRAY. The 



emperical thresholds may be printed out (plotted) if the proper 
option is input to the DISPLAY processor. 

When in the dot processing mode, however, the emperical threshold 
computations are not available to the user, due to the utilization 
of the training field storage areas for dot data storage, and due 
also to the fact that test fields are not allowed to be input 
when in the dot processing mode ( designated fields may be input, 
however) . 

This constitutes a revision to the output normally available via 
DSPLAY. When not in the dot processing mode (no DOTFILE control 
card input) , the usual emperical threshold options are available 
to the user. 

Also, when not in the dot processing mode, a call to PRTPCT is 
made from DSPLAY. This results in classification performance 
summaries by subclass, class, and category for either test 
fields or training fields. 

When in the dot processing mode, the test or training field 
performance summaries are not available to the user, due to the 
usage of the test and training field storage arrays described 
above . 

The test or training field performance summaries are replaced by 
dot classification performance summaries, when in the dot 
processing mode in the DISPLAY processor. 

F.7.5 STORAGE REQUIREMENTS 

DSPLAY requires 224 locations in the Insi ruction Bank and 40 
locations in the Data Bank for a total of 264 storage locations. 





F.7.6 DESCRIPTION 


DSPLAY is modified to bypass the call to the emperical threshold 
computation routine, EMTHRS, when in the dot processing mode 
(DOTKEY > 0 or DOTERR > 0) , and to bypass the call to the test/ 
training field classification performance output routine, PRTPCT, 
when in the dot processing mode. The modifications were made by 
inserting tests on the dot processing flags, DOTKEY and DOTERR, 
at the two appropriate places in DSPLAY. 

This modification is made necessary because the dot data utilizes 
the usual training field storage areas, and test fields are not 
allowed to be input to the DISPLAY processor when dot performance 
summaries are requested via the input DOTFILE control card. 

F . 7 . 7 FLOWCHART 
N/A 


F.7.8 LISTING 

The listing is available in the Data Techniques Laboratory, JSC 
Building 17. 


. 


F. 8 SOFTWARE SUBPROGRAM NO. 8 TAPLAB) 


TAPLAB is an assembly language (SLEUTHII) program for the Univac 
1108/EXEC 2, which has a Fortran V calling sequence so that it 
may be called from a Fortran V program. The purpose of TAPLAB 
is to provide the operational label (i.e. tape reel nvimber, or 
other label used on an "ASG" card in the equation "UNIT=LABEL") 
which is assigned to a logical tape unit (A,B,C,...). 

F.8.1 LINKAGES 

Subprogram TAPLAB presumes the availability of the UNIVAC EXEC 2 
system procedure (routine) , TLABL$. 


F . 8 . 2 INTERFACES 

Subprogram TAPLAB interfaces with the Fortran V calling program 
via two calling arguments. Subprogram TAPLAB references the 
UNI VAC /EXEC 2 system procedure, TLABL$. 


F.8.3 INPUTS 


Calling Sequence: CALL TAPLAB (UNIT, LABEL) 

Parameter Dimension Type In/Out 

UNIT 1 Integer In 


Definition 

The Fortran unit no. - 
legitimate numbers are 
(positive) 1-4,7-16, 
and 18-29. Fortran 
units 5,6 are system 
input, output units, 
respectively. Fortran 
unit 17 is reserved 
as system plot unit. 
Fortran unit 30 is 
the "reread" unit. 
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LABEL 1 


Alphanu- Out The operational 

meric or label assigned to 

INTEGER O'S UNIT. If UNIT=5,6, 

17, or 30 LABEL is 
set = XXXXXX and 
returned. If UNIT 
is unassigned as a 
tape unit, LABLE is 
returned all integer 
zero's in the six 
character positions 
in the word. If 
UNIT=l-4, 7-16, or 
18-29 and the 
corresponding logical 
unit (A-D, E-N, 0-Z) 
was assigned a tape 
unit by the ASG 
UNIT=LABEL card, 
LABEL = input tape 
label is returned 
as a 6-character 
alphanumeric . 


F.8.4 OUTPUT 

Subprogram TAPLAB provides to the calling program the tape label 
associated with (assigned to) the input (via first argument) 
Fortran unit number. The tape label is output via the second 
calling argument. The output tape label is either a six-character 
alphanumeric if the Fortran unit has been assigned a label, or 
all X's if the Fortran unit is unassigned, or all zero's if 
the Fortran unit is either 5,6,17, or 30. 



F.8.5 STORAGE REQUIREMENTS 

Subprogram TAPLAB requires 68 storage locations. 

F.8.6 DESCRIPTION 

TAPLAB testa the number in ARGl for within range 1 to 4. If 
YES, the unit number is converted to the appropriate logical 
unit character A - D , and a j\imp made to the call to the UNI VAC 
System Procedure, TLABL$. If NO, the number is tested for 
within range 7-16. If YES, the number is converted to the 
appropriate logical unit character E N and a jump made to 
the call to TLABL$. If NO, the number is tested for within 
range 18-30. If YES, the appropriate conversion to logical 
unit 0 - Z is made, and a jxmip to the call to TLABL$. If the 
unit is none of the number 1-4, 7-16, or 18-30, a jump is made 
directly to ERUNIT where a label = XXXXXX is placed in ARG2 
and a return made to the calling program. 

If the unit is a legitimate number within the ranges described 
above, but was not assigned to a tape unit in the run, by means 
of the UNIVAC ASG UNIT=LABEL control card, the return from 
TLABL$ is all zero's, and this is the label placed in ARG2 to 
return to the calling program. 

A test is made on the first character of the label returned to 
TAPLAB from TLABL$ . If the first character position contains 
a zero, a blank character is inserted in the first character 
position of the label returned via ARG2. This is for the case 
of a 5-character tape label having been assigned to the Fortran 
unit. 

TAPLAB utilizes registers AO, A7, A8, A9 , and Xll. 

F . 8 . 7 FLOWCHART 
N/A 





F.8.8 LISTING 


The listing is availabel in the EOD-LARSYS listing in the Data 
Techniques Laboratory, JSC Building 17. 
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G.l SOFTWARE SUBPROGRAM NO. 1 (TAPERD) 


The new capabilities of routine TAPERD are; 

1. on option, unpack sun agles from the header record of the 
MSS Universal format input tape. 

2. unpack the start and stop pixel number, pixel skip number, 
line skip number from the header record of the MSS Universal 
format input tape. 

G.1.1 LINKAGES 

The routine TAPERD is called by the routines; DOTS, STOMAP , 
RDDATA, CLSFY2, LEARN, NDHSTl , LNTRAN , HSTGRM, HISTGM. 

The corresponding processors are DOTDAT, LABEL, ISOCLS, CLASSFY, 
STAT, NDHIST, DATATR, HIST, and GRAYMP. 

G.l. 2 INTERFACES 

The subprogram TAPERD interfaces with other routines through 
common block ISOLNK. 

G.l. 3 INPUTS 

TAPERD has 3 entry points; TAPHDR, FLDINT and LINERD 
• Calling sequence; CALL TAPHDR (DATAPE, DATFIL) 


Parameter 

Dimension 

In/Out 

Definition 

DATAPE 

1 

In 

Tape unit number. 

DATFIL 

1 

In 

Number of EOF"s tc 
be read over in 
positioning tape, 
minus one. 

« Calling 

sequence ; CALL FLDINT 

(FLDINF, FETVEC, NOFEAT) 

Parameter 

Dimension 

In/Out 

Definition 

FLDINF 

6 

In 

Rectangular field 


description for use 




Par 2 irtieter 


Dimension 


In/Out 


Definition 


FLDINF 

(Continued) 


FETVEC 30 In 

NOFEAT 1 In 


• Calling sequence: CALL LINERD 

Parameter Dimension 

IDATA Depends upon 

buffer size 
in calling 
routine 

ENDTAP 1 Out 


G . 1 . 4 OUTPUTS 
N/A 

G . 1 . 5 STORAGE REQUIREMENTS 

Code - 16416g Data - 2420^ 

G.1.6 DESCRIPTION 

The information to be unpacked from the Universal 
occupies the following byte positions on the tape 
Byte Numbers Description 


(IDATA, ENDTAP) 
In/Out 

Out 


in unpacking 
pixels from 
MSS tape. 

Contains channels 
requested. 

Niomber of channels 
in FETVEC. 

Definition 

Array containing 
unpacked data. 

EOF trigger : 
returns value 
of -1 if an 
EOF is encountered 


format MSS tape 
header record: 


109-110 


pixel start number 


Byte Numbers 
110-111 
1789-1790 
1791-1792 
2201-2202 

Tape 2203-2304 
Sun Angles etc . 


Description 
pixel stop number 
pixel skip factor 
line skip factor 
sun angle for pass 1 
(channels 1-4) 
sun angle for pass 2 
(channels 5-8) etc. (up to 8 
passes ) 


This information is unpacked by simply extending the TAPHDR 
routine data vectors HWRD, BIT and NB, expanding the dimension of 
the ID vector, and expanding an existing unpacking code block. 

The information listed above will be placed in labeled common 
ISOLNK. 

Please refer to the document "Universal Format Read Routine for 
UNIVAC 1108 Programs", Technical Memorandum, Project Number 1651V, 
August, 1973 for a full description of utility routine TAPERD. 


G.1.7 FLOW CHART 
N/A 

G.1.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 




G. 2 SOFTWARE SUBPROGRAM NO. 2 (RDDOTS) 


Routine RDDOTS reads information from the dot data tape. 

G . 2 . 1 LINKAGES 

Routine RDDOTS is called by routines FILERD, ISOCLS, and SETUPS. 
The corresponding processors are LABEL, ISOCLS, and DISPLAY. 
RDDOTS calls a file positioning routine FSBSFL. 

G . 2 . 2 INTERFACES 

Interfaces with other routines through the common block GLOBAL. 
From this block is obtained the Fortran dot data tape unit number 
DOTUNT and the file number DOTFIL. The default unit number is 
8 (UN I VAC number F) . 

G . 2 . 3 INPUTS 

Calling sequence: RDDOTS (DOTS, DOTVEC, TOTDT3, TYPSWT, SIZES, 




TOTDT2, NOCAT, 
DOFEAT, FETVEC, 
TOTVRT, FLDSAV, 

CATNAM, NOFET2, FETVC2, 
NOSUN, ANGLE, NOFLD, 
VERTEX, KVAR) 

Parameter 

Dimension In/Out 

Definition 

DOTS 

(SIZES, 

TOTDT2) Out 

Dot information, 
spectral, spatial, or 
both depending upon 
option value of TYPSWT 

DOTVEC 

TOTDT 3 

In 

Depends upon valve of 


TYPSWT : 

1 - list of dots for 

which spectral 
information is 
wanted 

2 - ignored, except for 

a check if T0TDT3 

+ 0 . 




Parameter 


Dimension 


In/Out 


Definition 


DOTVEC (cont) 3 - list of dots to be 

excluded, all other dots 
on the DOTFIL are passed 
back to calling program 

T0TDT3 1 In Depends upon value of 

TYPSWT : 

1 - number of dots 

requested 

2 - should be preset to 

0 

3 - number of dots to be 

excluded 

TYPSWT 1 In Option switch, must be 

preset to a value of 
1,2, or 3 

1 - spectral information 

requested 

2 - spatial information 

requested 

(sample, line, type, 
and category numbers) 

3 - both spatial and 

spectral information 
requested 


SIZES 

1 

Out 

Depends 
TYPSWT ; 

upon value of 




1 - set 

to 

N0FET2 




2 - set 

to 

4 




3 - set 

to 

4 -i N0FET2 

T0TDT2 

1 

Out 

Depends 
TYPSWT ; 

upon value of 




1 ~ set 

to 

T0TDT3 




2 - set 

to 

TOTDOT 



Parameter 

Dimension 

In/Out 

Definition 

T0TDT2 

(continued) 

1 

Out 

(number of dots on 
DOTFIL) 

3 - set to TOTDOT-TOTDT3 

NOCAT 

1 

Out 

Number of categories, 
value extracted from 
DOTFIL 

CATNAM 

NOCAT 

Out 

Category names , values 
extracted from DOTFIL 

NOFET2 

1 

In/out 

Number of channels for 
which information is 
requested. If set to 0, 
reset to NOFEAT. If 
TYPSWT = 2, the value 
of NOFET2 is ignored. 

FETVC2 

NOFET2 

In 

Channel numbers for 
which information is 
requested, ignored if 
TYPSWT=2 . 

NOSUN 

1 

Out 

Number of sun angles 
on DOTFIL 

ANGLE 

NOSUN 

Out 

Sun angle values 
on DOTFIL 

NOFLD 

1 

Out 

Number of fields on 
DOTFIL 

TOTVRT 

1 

Out 

Number of vertices on 
DOTFIL 

FLDSAV 

(4, NOFLD) 

Out 

Field information on 
DOTFIL 


Q-6 '.';y7 


Parameter 


Dimension 


Definition 


In/Out 

(2, TOTVRT) Out Vertex information on 

DOTPIL 

(SIZES, TOTDT2) Out Set if TYPSWT = 1, 

same information as 
DOTS, but in floating 
point. 

G . 2 . 4 OUTPUT 
N/A 

G . 2 . 5 STORAGE REQUIREMENTS 
Code - 762g Data - 11767g 

G.2.6 DESCRIPTION 

The dot data tape is read and appropriate information is extracted. 
There are three modes of extraction: 

1. dot spectral information 

2 . dot spatial information 

3. both spectral and spatial information. 

In all cases, the count information (first record) is made avail- 
able. 

Each file on the dot data tape consists of three records. Please 
refer to the flow chart provided in Figure G.2 for details con- 
cerning contents of these records, also the file description in 
Appendix H. 

Dot spectral information consists of radiance values for those 


VERTEX 

KVAR 




dots and channels specified. Dot spatial infoirmation consists 
of : 


1. sample number 

2 . line number 

3. type number (1, for label/starting dots, 2, for bias dots) 

4 . category number 

for those dots specified by the user. 

G . 2 . 7 FLOW CHART 

See Figure G- 1 for flow chart of the RDDOTS routine. 

G. 2,8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 



RDOOTS 




Don, DOWEC, TOTtvrl, TYPSWT, 
SIZES, TOTDT2, HOCAT, CAltlAH, 
IHOfETZ, FfcTVCJ, MOFh;AT, ELTWr, 
MOSIW, M«CtE, MOri.0, TOTVBT, 
^rU)5AV, VERTEX, KYAR 

- .1 

( REWIKD DOTVNT ■'A 

CAU, rsssn, to v 

rOSITlON TAPE 


\ READ FIRST RECORD 

^NOCAT, NOFEAT, HOPLD, TOTVRT, TOTDOT, 

. NOSUN, CATNAM(I),I>1, NOCAT, SIZE 





C 


ji NO V ,/cALl. ERROR 
TOTDTJ ' 


RECORD TOTAL 
*T0TDOT_,TEMDOT 


<[TYPSW « I'N. 



Figur« G-f«* flowchart for Subproorami RDDOTS 


arc 


or/ginal page is 
POOR QUAUTY 







I FOR J » 1, SlRESl 
I i JJ » mVC3 (Jl 
j JJJ • (kk'l) * SHE + JJ -f 4 
. KVAR tJ.kl - FLOAT (TEMDOT {JJJl) 
j DOTS CJ,)c) - TEKDOT (JJJ) 

: SPECTRAL IMFOHATIOH 



'for k • 1, SISESi 

IdOTS Ik, II » TEMDOT (kk-fk) 

SPATIAL INFORMATION 

I 

i 



MAJOR LOOP: 

FOR k - 1, TOTDOTt 
IF THIS DOT IS EXCLUDED, BUMP kk ! 
, AND CO TO END ca- MAJOR LOOP. 


OTHEBWIFE, BUMP JJ, COMPUTE JJJ - 
i (k-1) *SIBE ( 

FOR J • 1, 4: j 

' DOTS (J,JJ) » TEMDOT (JJJ+J) ; 

FOR J-1, NCfTTJ: t 

I • FETVC3 {JIDOTS (4+J , JJ) »TEMDOT 
I (JJJtI+4) 




Flgura C't-- Continued 
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G.3 SOFTWARE SUBPROGRAM NO. 3 (WRTDOT) 


Subroutine WRTDOT writes the dot data tape. One file is written 
per call to WRTDOT. 

G.3.1 LINKAGES 

Routine WRTDOT is called by routine DOTS in the DOTDATA processor. 
WRTDOT calls the UNIVAC file positioning routine FSBSFL. 

G.3. 2 INTERFACES 

No common blocks are provided for WRTDOT. 

G.3. 3 INPUTS 

Calling sequence: WRTDOT (TOTDOT, NOSUN, FLDSAV, VERTEX, ANGLE, 

DOTS, NOCAT, CATNAM, SIZE, NOFET2 , 
TOTVT2, NOFLD2, UNIT, FILE) 


Parameter 

Dimension 

In/Out 

Definition 

TOTDOT 

1 

In 

Total number of dots to be 
written on DOTFIL 

NOSUN 

1 

In 

Number of sun angles to 
be written. 

FLDSAV 

(4,NOFLD2) 

In 

Field information to be 
written 

VERTEX 

(2, TOTVT2) 

In 

Vertices to be written 

ANGLE 

NOSUN 

In 

Sun angles to be written 

DOTS 

(SIZE, TOTDOT) 

In 

Dot information to be 
written 

NOCAT 

1 

In 

Number of categories 

CATNAM 

NOCAT 

In 

Category names 

SIZE 

1 

In 

4 + NOFET2 

NOFET2 

1 

In 

Number of features (channels) 
to be written on DOTFIL 




Parameter 

Dimension 

In/Out 

Description 

FETVC2 

N0FET2 

In 

Channel numbers 

T0TVT2 

1 

In 

Total number of 
vertices to be written 

N0FLD2 

1 

In 

Total number of fields 
to be written 

UNIT 

1 

In 

Unit number (Fortran) 
of DOTFIL 

FILE 

1 

In 

File number of DOTFIL 
for this file 


G . 3 . 4 OUTPUT 

An unformatted binary tape is written, with three records per file. 


G.3.5 

STORAGE 

REQUIREMENTS 

Code - 

562g 

Data 57g 

G.3.6 

DESCRIPTION 


The contents of the output tape files are described in Appendix 
H. All the information needed to write a file is brought in 
through the calling sequence to subprogram WRTDOT. 

G.5.7 FLOW CHART 

See Figure G-2 for flow chart of WRTDOT. 

G.3.8 LISTING 

Available in the Data Techniques Laboratory, JSC Building 17. 


G-12 


r 


1 


WRTDOT 

\ 


V 


^TOTOOT, HOeUN, FLDSAV, %RTf X, ' AMCCiT? 
DOTS, MOCAT, CATHAH, SIRE, NoJETZ, / 
rETVC2, TOTVT2, MCWL02, UMIT, TllM 


V 

/ \ 

REWIND DOTUHT^CAIE 
' rSBSn, TO POSITION 
' DOTUNT TAPE / 




Fiqurc C-2.- Flovrch.irt of Subprogram WRTDOT 




APPENDIX H 


DOT DATA FILE FORMAT 




iX)T DATA FILE FORMAT 
(DOTFIL) 


The tape, DOTFIL is output by the DOTDATA processor. The 
records are written with an unformatted FORTRAN write. 

A file is output for each TYPE of field(s). The file con- 
sists of the following records: 

field information 
data record 


TYPE 1 consists of labeling dots; TYPE 2 consists of bias 
correction dots. 

Rec. No. 1 




( Rec. No. 1 
\ Roc. No. 2 I 
Repeat for / 

each TYPE ^ Rec. No. 3 I 


E-O-F 


WRITE (UOTUNT) NOCAT, NOFEAT, NOFLD, TOTVRT, TOTDOT, NOSUN, 

(CATNAM(I), 1^1, NOCAT), SIZE 


Parameter 

Dimeiif. i on 

Def i nit ion 

NOCAT 

1 

Number of c itcvjory 



names . 

NOFEAT 

I 

Number of channels. 

NOFLD 

1 

Number of fields. 

'i'OTVR'i 

1 

Number of vertices. 

TOTDOT 

1 

Number of dots 

NOSUN 

1 

Number of sun angle 

CATNAM 

NOCAT 

Array containing th 



catecTory names. 

SIZE 

1 

4 + NOFEAT 


Rec. No. 2 


WRITE (DOTUNT) (FETVEC(I), I»l, NOFEAT) 

( (FLDSAVd, J) , I»l,4), J-1, NOFLD) , 

( (VERTEX (I, J) , J»l, TOTVRT) , 

(ANGLE (I), 1=1, NOSUN) 


Parameter 

Dimension 

Definition 


FETVEC 

NOFEAT 

Array containing 
channel numbers. 

the 

FLDSAV 

(4, NOFLD) 

Array containing 
field description 

the 

' • 

VERTEX 

(2, NOFLD) 

Array containing 
field vertices. 

the 

ANGLE 

NOSUN 

Array containing 
sun angles. 

the 


Rec. No. 3 


WRITE (DOTUNT) ( (DOTS (I, J) , 1=1, SIZE), J=l, TOTDOT) 


Parameter 


Dimension 


Definition 


DOTS 


(TOTDOT, SIZE) 



Array containing the 

dot information. 

DOTS ( 1 , 1 ) =sample number 
for dot X. 

DOTS ( 2 , 1 ) =-l ine nuj:iibor 
for dot i. 

DOTS (3, I) =type number 
for dot i. 

DOTS (4,1) =category 

number for 
dots i (optional) 


Parameter 


Dimension 


Definition 


DOTS (5,1) 

. =*dot 

. vector 

DOTS(4+NOFEAT,I) 
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